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Translocation  of  the  Global  Positioning  System  has  proved  to  be  a  highly 
accurate  method  of  position  determination  for  onshore  and  airborne  navi- 
gation, but  it  had  not  been  previously  evaluated  for  nearshore  hydro- 
graphic  surveys.  The  technique  of  translocation  for  hydrographic  opera- 
tions involves  the  simultaneous  reception  of  signals  from  the  GPS 
satellites  by  two  independent  receivers;  one  receiver  onboard  the  survey 
vessel  and  one  located  at  a  known  stationary  position.  A  position 
correction  is  obtained  at  the  stationary  receiver,  which  is  then  applied 
to  the  shipboard  receiver  online  or  during  postprocessing.  An  accuracy 
determination  of  translocated  GPS  was  conducted  at  the  Naval  Postgraduate 
School,  Monterey,  CA,  May,  1981.  Two  methods  of  positioning  were  used 
for  comparison  with  GPS:  1)  a  least  squares  solution  of  three  lines  of 
position  observed  from  three  Wild  T-2  theodolites,  and  2)  a  position 
determined  from  ARTEMIS,  a  range-azimuth  short  range  microwave  position- 
ing system.  Translocated  GPS  accuracies  of  10  meters  were  determined. 
It  is  anticipated  that  greater  accuracies  will  be  obtained  by  using  a 
more  sophisticated  receiver  and  more  advanced  processing  methods. 
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I.   INTROPgCTION 

Shipboard  operations  for  the  purpose  of  oceanographic 
research,  commercial  transport,  tactical  operations,  or 
hydrographic  surveying  require  accurate  and  dependable 
positioning  systems.  The  operation  of  most  electronic 
positioning  systems  is  hampered  by  atmospheric  effects  on 
the  propagated  signal  and  design  limitations  which  restrict 
operating  ranges,  positional  accuracy  and  dependability.  A 
universal  positioning  system  that  provides  accurate 
information  to  the  user  regardless  of  atmospheric  conditions 
and  that  is  free  of  range  limitations  is  desirable.  with 
the  advent  of  satellite  positioning  techniques,  such  a 
system  is  now  available. 

The  Navy  Navigation  Satellite  System  (NAVSAT) ,  or 
TRANSIT,  has  been  operational  since  1964  and  has  proved 
dependable  and  sufficiently  accurate  for  certain  users. 
Although  developed  for  military  uses,  the  application  of 
TRANSIT  in  the  civilian  community  has  continued  to  expand  at 
an  exponential  rate.  The  system,  however,  is  limited  by  the 
time  interval  between  fixes,  which  is  90  minutes  on  the 
average,  and  the  requirement  for  precise  Icnowledge  of  the 
ship's   velocity.    Although   the   system   does  not   provide 


16 


One  method  of  improving  system  accuracy  is  by  using  two 
receivers  simultaneously  in  a  technique  termed 
translocation,  or  differential  mode.  The  method  is  based  on 
the  assumption  that  signal  propagation  errors  to  two 
proximate  receivers  are  nearly  the  same.  One  receiver  is 
placed  over  a  known  position,  and  the  variation  of  the 
observed  position  with  time  is  determined  relative  to  the 
absolute  position.  This  variation  is  applied  to  the 
observed  positions  determined  by  a  simultaneously  operating 
receiver  located  at  another  location — on  a  survey  vessel, 
for  example.  The  assumption  underlying  this  technique  is 
that  the  distance  between  the  two  units  is  not  sufficient  to 
introduce  a  significant  additional  propagation  error.  The 
resulting  position  of  the  translocated  receiver  should  be 
more  accurate  than  the  position  obtained  using  a  single 
receiver. 

For  hydrographic  operations,  the  translocated  Global 
Positioning  System  (GPS)  provides  an  alternative  to 
presently  used  positioning  systems.  The  satellites  orbit  at 
20,000  1cm,  which  theoretically  would  allow  receiver 
separation  of  as  much  as  500  Icm  without  causing  significant 
additional  system   error.     With  this   flexibility,    one 
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sufficient  positional  accuracy  for  hydrographic  surveying, 
it  has  been  used  for  geodetic  application  by  the  Defense 
Mapping  Agency  (DMA)  using  fixed  shore-based  receivers  to 
establish  positioning  networks  in  inaccessible  areas 
worldwide.   [ Ref .  1] 

The  Global  Positioning  System,  or  NAVSTAR,  is  a  more 
sophisticated  system  designed  to  replace  TRANSIT  in  the 
mid-1980s.  Six  satellites  are  currently  in  orbit  with 
twelve  more  to  be  launched  before  the  operational  date  of 
1987.  The  18-satellite  orbit  configuration  will  provide  a 
multi-user,  passive  system  with  2t*-hour  availability,  global 
coverage,  and  possible  fix  intervals  of  less  than  one 
second.  The  operational  network  consists  of  three 
elements — the  satellites,  the  ground-based  control  segment, 
operated  by  the  Space  and  Missile  Systems  Organization 
(SAMSO) ,  and  the  user's  receiver.  Prototype  receivers  have 
been  tested  under  optimum  conditions  by  SAMSO  at  the  Yuma 
Proving  Ground,  Arizona,  on  both  mobile  and  stationary 
platforms.  Results  show  excellent  system  stability  and 
accuracy  when  compared  to  currently  used  positioning 
systems.   C^ef.  2] 
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shorebased  receiver  could  be  placed  at  a  central  location 
and  a  survey  of  several  hundred  kilometers  of  coastline  or 
offshore  area  could  be  completed  with  only  one  additional 
receiver  on  the  survey  vessel.  [Ref,  3]  The  value  of  this 
flexibility  would  be  realized  in  the  significant  decrease  in 
survey  planning,  cost,  and  time  compared  to  that  spent 
presently  in  conducting  hydrographic  operations. 

The  accuracy  of  the  system  would  be  sufficient  for  all 
survey  scales.  The  National  Ocean  Survey  (NOS)  defines  the 
required  accuracy  of  a  positioning  system  to  be  no  more  than 
one-third  of  the  survey  accuracy,  which  is  1.5  millimeters 
at  the  scale  of  the  survey.  [Sef.  4]  Non-translocated  GPS 
has  been  shown  to  be  of  sufficient  accuracy  for  1:80,000 
scale  surveys  and  smaller.   [Bef.  5] 

The  object  of  our  research  is  to  determine  whether  GPS 
at  the  present  operational  level  is  adequate  for  1:10,000 
scale  surveys.  To  do  this,  a  position  determined  by  two  GPS 
receivers  in  the  translocated  mode  is  compared  to  the  same 
position  determined  by  two  other  methods.  If  translocated 
GPS  can  be  shown  to  achieve  specified  inshore  survey 
accuracies,  the  system  would  provide  sufficient  accuracies 
at  all  smaller  scale  surveys. 
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Part  of  the  difficulty  in  ietermining  this  level  of 
accuracy  concerns  the  reference  system  used.  Most  geodetic 
positions  published  by  the  I?ational  Geodetic  Survey  (!IGS) 
are  based  on  the  North  American  Reference  Datum  of  1927 
(IIAD-27)         and        the      Clarice         1866      Ellipsoid,  the      best 

approximation  to  the  shape  of  the  earth  for  the  continental 
Onited  states.  The  satellite  systems,  however,  use  a 
reference  system  which  is  not  based  on  a  portion  of  the 
earth's   surface.  A   mass-centered   ellipsoid      referenced   to 

the  earth's  center  of  mass  has  been  computed  using  observed 
gravity  and  astronomical  data.  The  best  fit  to  this  global 
ellipsoid  is  the  World  Geodetic  System  of  1972  (MGS-72) , 
which  is  utilized  by  both  '' TRANSIT  and  GPS.  Transformation 
equations  exist  which  convert  NAD-27  coordinates  to  WGS-72 
coordinates,  but  uncertainty  in  the  gravity  model  produces 
errors  as  large  as  10  meters  due  to  the  transformation 
process   alone,      [Sef.    6] 

To  establish  all  positions  on  the  same  reference 
ellipsoid  for  this  research,  a  TRANSIT  receiver  was  acquired 
from  DMA  and  a  first  order  Doppler  station  was  established 
at  the  central  location  of  the  test  area.  A  third  order 
geodetic   survey     was   conducted    from     the    TRANSIT      station   to 
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establish  two  additional  stations  on  the  WGS-72  ellipsoid. 
Two  single  channel  Manpack  GPS  receivers  were  obtained  from 
SAMSO.  Stability  testing  and  a  receiver  comparison  test 
were  conducted  prior  to  shipboard  testing  to  determine  a 
guanxitative  measure  of  receiver  performance.  0ns  of  the 
receivers  was  then  placed  over  the  TRANSIT  station  mark  and 
the  other  aboard  a  survey  vessel  which  maneuvered  along 
predetermined  track  lines.  Data  were  recorded  via  an 
interface  to  an  HP-9825  computer  and  printer  at  each  site. 
Simultaneously,  two  additional  independent  measurements  were 
made  of  ship  position — one  using  a  range-azimuth  short 
range  microwave  positioning  system  (ARTEMIS)  and  the  other  a 
position  determined  by  a  least  squares  solution  of  three 
lines  of  position  observed  from  three  1-second  theodolites 
set  on  the  three  pre-established  WGS-72  locations.  The 
comparison  of  the  ARTEMIS  and  theodolite-derived  ship 
positions  with  the  translocated  GPS  ship  position  provided 
the  test  results. 
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II.       EOaiPMENT    DESCRIPTION    AND    OPERATIONAL    CAPABILITIES 
Two   positioning    systems   were     chosen   for   comparison   with 
the  Global   Positioning      System:      a   theodolite    network      and  a 
microwave      positioning      system.  The     theodolite      network 

consisted  of  three  shore-based  units  so  as  to  provide  three 
lines     of      position.  The        microwave      positioning      system 

provided, a  range  and  an  azimuth  from  a  known  location  on 
shore.  Both   systems      are     accurate     at   short     ranges     and 

represent  two  contrasting  methods  of  positioning--visual  and 
electronic.  The  theodolite  method  of  positioning  at  short 
range  has  been  and  continues  to  be  used  for  hydrographic 
operations  as  an  alternative  to  electronic  systems  when 
geometry  and  system  accuracy  limitations  impose  the  need  for 
a  more  flexible  positioning  system.  The  requirement  for 
experienced        theodolite        observers        as  well        as        the 

communication  and  logistics  difficulties  involved  limit  the 
use   of   this  method   to   special    cases. 

This  chapter  describes  the  ARTEMIS  short-range 
electronic  positioning  system,  the  operation  of  the 
theodolite  network  as  used  during  the  comparison  test,  and 
the  TRANSIT   satellite  system. 
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i,       ARTEMIS  POSITIONING  SYSTEM 

The  ARTEMIS  position  fixing  system  is  a  short-range 
licrowave  system  (to  30  km),  bailt  by  Christiaan 
luygenslaboratorium  B.  V.,  in  the  Netherlands,  It  has  been 
ised  operationally  (primarily  in  Europe)  since  1972,  when 
:he  system  was  introduced  at  the  International  Hydrographic 
:onference  in  Monaco.  At  the  present  time,  fifty  systems 
ire  in  use.  The  system  is  unique  in  design,  is  highly 
iccurate  at  short  ranges,  and  is  particularly  useful  for 
small  harbor  surveys,  river  surveys  and  relative  positioning 
ises.  For  this  study,  the  ARTEMIS  system,  along  with 
technical  assistance,  was  obtained  from  the  American 
iistributor,  MARINA7  Corporation  of  Houston,  TEXAS. 

ARTEMIS  is  a  range-azimuth  positioning  system  which 
employs  a  unique  tracking  method.  A  stationary  or  fixed 
unit  is  placed  over  a  known  position  on  shore.  Osing  a 
telescope  that  is  mounted  on  the  antenna,  the  antenna  is 
sighted  onto  another  known  point  and  the  known  azimuth 
between  the  two  locations  is  entered  into  the  system  via  a 
digital  display  on  the  fixed  unit.  The  unit  is  then  locked 
into  place  for  operation.  A  mobile  unit  is  placed  on  the 
vessel  whose  position  is  desired. 
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The  mobile  and  fixed  antennas  resemble  radar  antennas 
which,  while  operating,  track  each  other  so  that  the  two 
antennas  are  always  parallel  to  each  other  and  perpendicular 
to  the  line  of  position  between  them.  A  microwave  link  is 
established  between  the  two  antennas  during  periods  of 
measurement  and  information  is  relayed  between  the  mobile 
and  fixed  units.  (A  voice  link  is  also  available  during 
operation.)  The  fixed  unit  provides  the  direction  and  the 
azimuth  of  the  mobile  unit  with  respect  to  the  point  sighted 
on  by  the  telescope.  At  the  same  time,  the  mobile  unit 
transmits  a  coded  interruption  of  the  signal  to  the  fixed 
unit,  which  replies  with  the  same  interruption.  The  time 
difference  between  the  transmission  and  reception  at  the 
mobile  unit  determines  the  distance  between  the  two 
antennas. 

Distance  accuracy  is  determined  by  the  number  of 
measurements  averaged  over  time  and  is  available  in  two 
modes — a  dynamic  mode  which  averages  1000  measurements  or  a 
static  mode  which  measures  10,000  measurements  per  displayed 
distance.  The  azimuth  and  distance  information  is  displayed 
by  the  mobile  unit  so  that  personnel  are  not  required  at  the 
stationary  unit  once  operation  has  begun.   The  system  can  be 
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operated  in  long  or  short-range   nodes,   which  determine  the 
amount  of  power  output  for  the  system. 

The  microwave  freguencies  used  are  9.2  and  9.3  GHz. 
Problems  arising  with  signal  interference,  observed  with  all 
microwave  positioning  systems,  therefore  can  occur  with 
ARTEMIS.  The  regions  of  range  holes,  or  multipath 
interference,  can  be  computed  for  up  to  10  !cm  range  by  prior 
knowledge  of  the  antenna  elevations  above  sea  level  and  the 
operating  distance  from  the  shore- based  antenna.   [ Sef .  7] 

Sach  antenna  is  designed  with  a  22°  vertical  beam  width 
and  a  2°  horizontal  beam  width.  A  power  source  of  22  volts 
DC  per  unit  is  required,  with-  an  average  current  load  of  2.5 
amps. 

The  accuracy  of  the  system   depends  on  both  the  distance 
and  the   angle  measured.    The  angle   measurement  has  a  2  a 
probability  of  measuring  within  2'   of  arc   (0.033°) .    The 
distance   measurement  has  a  2  cr  probability   that  the   mean 
distance  value  will  vary  within  +  1.5  meters. 

One  of  the  advantages  of  ARTEMIS  over  other  electronic 
systems  is  that  the  two  lines  of  position  defining  the  fix 
always  intersect  at  90*^  .  A  unique  advantage  is  that  the 
coverage   around  the   fixed  unit   can  be   circular,   with  a 
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aaximum  range  cf  operation  to  30  km,  dependent  on  antenna 
heights.  The  short  range  applications  of  the  ARTEMIS  fill  a 
need  not  provided  by  many  currently  used  positioning 
systems. 

For  this  study,  the  ARTSMIS  was  used  in  conjunction  with 
the  positions  derived  from  the  theodolite  network  to 
determine  the  best  estimate  of  the  ship's  position.  The 
accuracy  of  the  ARTEMIS  was  more  than  adequate  for  the 
ranges  encountered  during  the  test,  and  the  ease  of  use  and 
operational  versatility  of  the  system  made  its  use  very 
advantageous. 

B.   THEODOLITE  NETWORK 

A  network  of  three  Wild  r-2  theodolites  was  used  as  the 
control  positioning  system.  The  theodolite  is  widely  used 
for  surveying  as  a  means  of  measuring  angles  between  two 
observational  locations — either  in  a  horizontal  plane  for 
traverse  or  triangulation  surveys,  or  a  vertical  plane  for 
astronomical  observations. 

The  use  of  three  theodolites  as  a  short-range 
positioning  system  was  an  accepted  method  of  positioning 
prior  to  development  of  electronic  distance  measuring 
systems.    There  remain  disadvantages   when  using  electronic 
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positioning  systems  for  hydrogra  phic  operations  at  short 
ranges  (1-2  km) .  Most  short-range  electronic  systems  are 
designed  for  measurement  of  position  at  line-of-sight 
ranges.  This  region  includes  a  large  portion  of  inshore 
hydrographic  surveys.  Inshore  positioning  involves  more 
thorough  planning  to  obtain  the  required  accuracy  for  the 
scale  of  the  survey.  For  example,  in  a  harbor  or  cove  it 
may  be  difficult  to  place  two  range-range  electronic  units 
so  that  the  area  to  be  surveyed  is  not  within  the  region  of 
positioning  system  uncertainty,  or  30°  (radially)  on  each 
side  of  the  baseline  between  the  two  units.  To  reposition 
the  electronic  units  as  often  as  is  required  to  obtain 
optimum  geometry  is  time-consuming  and  requires  manpower  and 
extensive  planning,  and  can  result  in  operational  delays. 
Interference  of  the  signal  due  to  reflection  from  a  calm  sea 
surface,  or  multipathing,  is  also  a  problem  at  short  ranges 
for  microwave  systems. 

Theodolites  can  be  set  at  any  intervisible  location  with 
preplanned  geometry  and  can  be  used  for  the  entire  region  to 
be  surveyed.  This  method  requires  less  expensive  equipment 
(three  theodolites) ,  more  manpower  (a  recorder  and  an 
observer   at  each  site) ,   communication   between  the   three 
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locations  (via  portable  radio,  for  example)  and, reasonable 
w€ather  conditions.  The  basis  of  this  procedure  is  that 
three  angles  observed  from  knovn  locations  to  a 
predeterained  point  on  the  vessel  constitute  three  lines  of 
position.  After  removing  all  possible  systenatic  errors  and 
through  a  least  squares  solution  of  the  three  lines  of 
position,  the  position  of  the  vessel  is  determined.  A  major 
disadvantage  is  that  observational  errors  are  not  determined 
in  real  time  because  the  data  cannot  be  processed 
imnediately. 

The  basic  procedure  for  the  method  requires 
predetermined  sites  which  are  known  relative  to  each  other. 
This  may  be  achieved  by  using  published  geographic  positions 
or  by  conducting  a  local  plane  survey  to  establish  new 
positions  using  preestablished  stations  nearby.  Either  way, 
all  three  locations  must  be  based  on  the  same  datum  and 
preferably  permanently  positioned  by  a  disk  or  a  station 
mark. 

Each  theodolite  was  aligned  using  a  plumb  bob  or  optical 
plummet  over  one  of  the  known  station  marks  or  disks.  3y  a 
predetermined  scheme,  each  site  was  marked  so  as  to  be 
clearly  visible  to  the   other  sites   (lights  were   used  for 
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this  study) .  The  telescope  of  each  theodolite  was  sighted 
onto  one  of  the  other  sites  and  the  horizontal  circle  set  to 
read  an  arbitrary  value  close  to  O''.  The  crosshairs  in  the 
telescope  were  set  to  intersect  the  aark.  When  the 
theodolite  was  exactly  sighted,  the  horizontal  circle  was 
read  and  recorded.  (At  this  time  the  lights  were  turned 
off.)  The  operation  began  when  all  three  sites  completed 
this  step,  called  the  "initialling  procedure",  or  "initial 
pointing".  The  instrument  was  sighted  on  a  point  on  the 
survey  vessel  that  was  visible  to  all  operators.  The  ship 
was  trac)ced  by  the  operators  by  using  the  micrometer  wheels 
on  the  instrument  until  a  "mark"  was  called  over  the  radio 
by  the  recorder  at  the  control  theodolite  location.  At  that 
instant,  the  operator  stopped  tracking  and  the  horizontal 
circle  was  read  and  recorded.  This  observation  constituted 
a  fix.  After  a  series  of  fixes  wers  recorded,  or  when  the 
survey  vessel  reached  a  predetermined  location,  a  set  of 
fixes  was  completed.  The  lights  were  then  reactivated  at 
each  site,  each  observer  intersected  the  original  mark,  the 
horizontal  circle  was  again  read  and  the  "final  pointing" 
was  recorded.  If  necessary,  the  instrument  was  relevelled. 
The  process  was  then  repeated  by  starting  with  the 
initialling  procedure. 
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It  was  assumed  that  the  theodolites  were  in  good 
operating  condition.  Collimation  error  had  previously  been 
determined.  The  observers  were  kept  at  the  same  locations 
throughout  the  test  to  minimize  the  error  due  to  individual 
observer  variations.  Since  only  one  portion  of  the 
horizontal  circle  was  used,  and  no  reverse  readings 
observed,  a  systematic  error  was  introduced.  This  was 
assumed  to  be  minimized  during  computation  of  the  vessel 
position. 

The  observations  were  corrected  for  collimation  error 
and  for  pointing  error  by  evenly  distributing  the  difference 
between  the  initial  and  final  pointing  values.  The 
resulting  angles  were  processed  using  an  intersection 
computation  with  a  least  squares  fit  applied  to  optimize  the 
result.  The  computations  for  this  study  were  processed 
using  an  HP-9815  computer  and  software  prepared  by  Cdr 
Ludvig  Pfeifer,  NOAA.  The  resulting  position  was  used  as 
the  best  estimate  of  the  position  of  the  ship. 

C.   THANSIT  POSITIONING  SYSTEM 

The  TRANSIT,  or  Navy  Navigation  Satellite  System,  is  the 
most  widely  utilized  satellite  positioning  system  at  the 
present  time.   The  network  of  satellites  and  ground  tracking 
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stations  has  been  in  operation  since  1964  and  is  becoming 
widely  used  by  both  the  military  and  civilian  communities. 
It  is  expected  to  continue  operation  until  at  Isast  1995, 
when  the  Global  Positioning  System  is  expected  to  replace 
it. 

The  system  consists  of  a  series  of  TRANSIT  satellites 
(fig-  2.1.)  in  fixed  polar  orbits  of  1095  km  elevation, 
forming  a  network  within  which  the  earth  rotates  (fig. 2. 2.). 
The  orbits  are  spaced  so  that  satellite  passes  occur,  for  a 
given  ground  location,  at  intervals  of  between  35  and  100 
minutes.  Each  satellite  travels  an  observable  distance  of 
between  4400  and  7000  km  per  pass,  providing  a  sufficient 
time  interval  and  baseline  over  which  the  user  can  record 
data. 

The  idea  of  positioning  an  observer  on  the  earth's 
surface  by  using  the  satellite  as  a  reference  system 
originated  from  observation  of  Sputnik  1,  launched  by  the 
Soviet  Onion  in  1957.  By  recording  the  change  in  frequency 
of  the  satellite  signal  as  the  satellite  passed  overhead,  a 
Doppler  shift  could  be  measured.  The  Doppler  shift  gives  a 
spatial  relationship  as  a  function  of  time  between  the 
observer   and   the   satellite.     If  the   position   of   the 
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Fig.   2.1.      TRANSIT  SataKite 
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"ig.  2.2.  Transit  Satellite  Orbital  Network 
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satellite  is  known  at  the  tiae  of  observation,  the 
observer's  position  can  be  determined. 

In  the  TRANSIT  system,  the  satellite  positions  are 
precisely  determined  by  tracking  stations  located  at 
Prospect  Harbor,  Maine;  Rossnount,  Minnesota;  and  Wahiawa, 
Hawaii.  The  tracking  information  (in  the  form  of  the 
Doppler  frequency  shift  as  a  function  of  time)  is  relayed  to 
the  computing  center  at  the  Naval  Astronautics  Group 
Headquarters,  Pt.  Mugu,  California.  Here  the  actual 
satellite  orbits  are  computed.  This  information  is  used  to 
predict  orbit  information.  Parameters  of  this  predicted 
orbit  are  relayed  to  two  injection  stations  at  Rosemount  and 
Pt.  Mugu.  During  the  next  satellite  pass,  the  information  is 
relayed  to  the  satellite  as  part  of  the  navigation  message, 
with  new  updates  occurring  every  12  hours. 

The  satellites  therefore  contain  predicted  orbit 
information  which  is  received  by  the  user  via  a  satellite 
receiver.  Several  receivers  are  available  with  complexity 
and  positional  accuracy  dependent  upon  the  needs  of  the 
user,  A  computer  is  normally  required  to  determine  the 
user's  mobile  or  fixed  position  in  real  time  from  a  single 
satellite  pass. 
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For  many  users  having  a  dual  channel  receiver  and  known 
velocity,  this  method  of  positioning  provides  sufficient 
accuracy  with  typical  error  of  27  to  37  meters  root  sum 
square.  Optimum  accuracy  is  obtained,  however,  if  the 
receiver  position  is  fixed,  and  several  consecutive  passes 
of  the  same  satellite  are  recorded.  During  each  pass, 
lasting  from  10  to  16  minutes  as  the  satellite  travels  from 
one  horizon  to  the  other,  many  Doppler  counts  (observations 
of  shift  in  frequency  as  a  function  of  time)  are  recorded. 
After  a  predetermined  number  of  passes  are  stored,  the  data 
can  be  forwarded  to  the  Defense  Mapping  Agency  in  Brookmont, 
Maryland.  There  the  precise  ephemeris  data  from  the  ground 
tracking  system  is  applied  to  the  data  from  the  recorded 
passes  to  provide  the  optimum  solution.  The  resulting 
position  has  a  typical  error  of  6.3  meters  root  mean  square 
for  a  series  of  passes,  with  a  3-dimensional  result  of  1.5 
meters  per  axis  repeatability  after  25  precise  ephemeris 
passes.  This  method  is  used  to  determine  positions  in 
r smote  areas  by  DMA  geodetic  survey  groups  using  the 
Magnavox  AN/PSR-14  Geoceiver  (fig. 2. 3.).   [Hef.  8] 

The  TRANSIT  system  uses  the  World  Geodetic  System  of 
1972  (WGS-72)    reference  ellipsoid   which  is   based  on   the 


35 


^^ 


h 


X 


"7 


^ 


i 


AN/PRR-14  GEOCEIVER 
ANTENNA 

Fig  .  2.3  . 


35 


Universal  Space  Rectangular  Coordinate  System — an  earth, 
mass-centered  coordinate  system.  The  WGS-72  ellipsoid 
describes  the  best  fit  of  the  entire  earth  to  a  mathematical 
model.  The  satellite  orbits  are  computed  using  this 
geopotential  model,  the  accuracy  of  which  is  based  on 
gravitational  and  astronomical  data. 

For  the  present  study,  a  Magnavox  AN/PRR-14  Geoceiver 
was  acquired  from  the  Department  of  Satellite  Geophysics, 
Satellite  Tracking  Branch,  DMA.  The  receiver  was  placed 
over  a  fixed  point  and  t*0  passes  of  one  satellite  were 
observed.  The  procedures  described  previously  for  an 
optimum  solution  of  point  positioning  were  followed.  The 
resulting  latitude,  longitude,  and  height  relative  to  the 
WGS-72  ellipsoid  are  coordinates  of  a  first  order  TRANSIT 
position  as  defined  by  DMA  standards.  This  position  was 
used  as  the  basis  for  horizontal  control  for  the  study  (see 
test  procedures) . 


37 


III.   GPS  SYST5I1  DESCRIPTION  AND  OPERATIONAL  CAPABILITIES 

A.   SYSTEM  DESCRIPTION 

NA7STAE,  or  the  Global  Positioning  System,  was  designed 
for  maximum  operational  capability — high  accuracy, 
twenty-four  hour  access,  passive  user,  and  world-wide 
coverage.  A  brief  description  of  the  GPS,  and  the 
satellites  (NAVSTARs),  will  serve  as  a  background  for  the 
more  in-depth  discussion  immediately  following. 

GPS  navigation  is  based  on  the  measurement  of  four 
ranges,  or  pseudo-ranges,  from  each  of  four  satellites  with 
known,  or  predicted,  positions.  Each  range  determination  is 
computed  using  the  propagation  velocity  of  the  signal 
multiplied  by  the  satellite-to-user  travel  time.  Satellite 
transmission  time  is  included  as  part  of  the  navigation 
signal  message.  Prom  the  four  known  ranges,  the  user 
position  and  user  clock  error  are  computed  on  a  three-axis, 
orthogonal  coordinate  system.  This  computation  is  then  used 
to  transform  the  user  position  onto  the  WGS-72  ellipsoid. 

The  remainder  of  this  chapter  will  provide  the 
background  for  GPS  development  and  the  three  :najor  system 
segments  of  the  GPS:   the  space  system  segment,   the  control 
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segaent,  and  the  user  system  segment.  Included  in  the  user 
system  segment  will  be  a  discussion  of  the  Jlanpack  receiver 
as  developed  by  Texas  Instruments  Corporation  used  for  this 
study. 

Hith  the  successful  operation  of  the  TRANSIT  system  in 
the  early  1960»s,  the  Department  of  Defense  (DOD)  began  to 
consider  the  specifications  of,  and  the  technology  required 
for,  the  development  and  implementation  of  a  second 
generation  system.  On  17  April,  1973,  the  Deputy  Secretary 
of  Defense  stated  that  a  Defense  Navigation  Satellite  System 
was  to  be  developed.  [Ref.  9]  Osing  the  new  synthetic 
oscillator  technology,  and  the  Department  of  Defense  surface 
and  nearsurface  navigation  objectives  and  requirements,  a 
DOD  interdepartmental  task  force  developed  the  GPS  as  the 
second  generation  satellite  navigation  system.  With  the 
United  States  Air  Force  (USAF)  as  the  managing  agency,  the 
GPS  became  a  joint  service  program,  with  representatives 
from  the  DMA,  Army,  Navy,  Marines,  Coast  Guard,  and  the 
North  Atlantic  Treaty  Organization-  [Ref.  10]  The  Joint 
Programs  Office  was  established  within  the  Air  Force  Space 
and  Missile  Systems  Organization  to  design,  develop,  and 
iipiement  the  GPS  program. 
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GPS  was  conceived  and  designed  to  provide  continuous, 
real-time,  world-wide  navigation  coverage  for  subsurface, 
surface  and  nearsurface  operational  vessels,  vehicles,  and 
aircraft.  Accuracy  of  better  than  10  meters  was  desired 
under  all  conditions,  with  navigational  capabilities  to 
update  and  interface  with  other  forms  of  existing  and 
projected  navigational  systems.  Another  requirement  of  the 
GPS  satellites  was  nearly  trouble-free  and  simple  user 
operation,  even  in  low  signal-to-noise  environments.  An 
encoded  signal  structure  was  required  for  the  satellite 
system  for  security  reasons.  In  addition,  the  signal 
transmitted  by  the  NA7STAR  satellites  had  to  include  other 
operational  information  in  conjunction  with  the  navigational 
message.  To  provide  for  navigational  requirements  in  event 
of  hostilities,  the  GPS  was  designed  to  permit  the  use  of 
redundant  satellites. 

To  achieve  these  objectives  and  requirements,  the  GPS 
program  was  organized  according  to  three  developmental 
tasics:  concept  validation,  full-scale  engineering 
development  and  systems  tests,  and  production  and  operation. 
Starting  in  1974,  the  concept  validation  task  was  initiated 
with  the  merger  of  the  TIMATION  series  satellite  program  and 
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the  OSAF  system  21B  (a  highly  accurate  three-diiensional 
navigation      system)       program.  Tests      with   the     Navigation 

Technology    Satellite,      NTS-1,      utilizing    two   rubidium    clocks 

1  3 

(accuracy  5-10  parts  per  10  )  and  an  orbit  7500  nautical 
miles  (nm)  above  the  surface  and  three  ground  based 
receivers  (called  inverted  range)  provided  concept 
verification.  [Ref,  11]  Additional  tests  were  conducted 
with  an  improved  NTS-2  satellite  positioned  at  an  altitude 
of  10,980  nm  and  containing  two  cesium  clocJcs  (accuracies  on 
the  order  of  1  part  per  10^^).  [Ref.  12]  Orbital  height  of 
10,980  nm  permitted  the  evaluation  of  the  effect  of  the 
planned  GPS  constellation  altitude  on  navigational 
accuracies  using  the  improved  cesium  frequency  standards. 
In  1980,  NTS-3  was  launched  and  the  Advanced  Developmental 
Model  (ADM)  hydrogen  maser  frequency  standard  evaluated  for 
use  on  the  GPS  (providing  an  expected  frequency  accuracy  of 
1  part  per  10^^).  [Bef.  13]  The  results  obtained  from  the 
NTS-3  satellite  have  not  yet  been  released.  Since  1978, 
additional  satellites  have  been  orbited  to  provide  the 
six-satellite  constellation  required  to  fully  evaluate  the 
first  development  task — concept  validation.  Phase  I. 
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During  Phase  1,  user  requirements  were  consolidated  to 
allow  the  development  and  construction  of  four  basic  types 
of  test  receivers:  continuous  traclcing  receivers, 
sequential  traclcing  receivers,  lew  cost  receivers,  and  the 
Manpaclc  receivers.  [  Ref .  11]  These  receiver  types  were 
then  constructed  for  Phase  I  testing  by  various  contractors 
and  their  performances  evaluated  aboard  helicopters, 
aircraft,  ships,  vehicles  and  personnel  at  the  Joint 
Programs  Office's  testing  facilities.  Phase  I,  concept 
validation,  was  completed  in  late  1979  and  Phase  II  testing 
was  begun  in  early  1980  with  full  scale  engineering, 
developmental,  and  systems  tests. 

Phase  II  tasks,  finishing  the  operational  and  control 
segments  of  the  testing,  are  scheduled  for  completion  in 
mid-1983.  The  operational  control  station  will  be  relocated 
at  a  central  continental  anited  States  site  during  this 
phase.  Additionally,  all  satellite  monitor  receivers  at  the 
remote  monitoring  stations  will  be  updated  to  operational 
status.  Oser  equipment  and  receiver  selection  will  be  made 
in  two  increments.  Four  contractors  will  be  asked  to 
compete  for  the  second  increment  of  the  receiver  selection 
process  by  providing  upgraded  and  refined  user  receivers  for 
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extensive  testing  and  evaluation.  Two  of  the  four 
contractors  will  then  be  selected  to  continue  the 
coipetition  and  will  provide  additional  receiver  design, 
refinement  and  modification,  leading  to  prototype  production 
and  extensive  field  testing. 

When  Phase  II  is  completed  in  mid-1983.  Phase  III, 
production,  will  commence  with  the  selection  of  one  of  the 
two  contractors  to  produce  the  user  receivers.  Manufacture 
of  these  receivers  will  commence  shortly  after  selection, 
and  commercial  availability  of  receivers  is  expected  in 
1985.  [Ref.  15]  The  GPS  should  be  fully  operational  with 
world-wide  coverage  using  a  constellation  containing  a  total 
of  eighteen  satellites  in  1988.   [Ref.  16,  173 

The  GPS  is  considered  to  be  composed  of  three 
inter-related  system  segments:  the  space  system  segment, 
the  control  segment,  and  the  user  system  segment.  Of  these, 
the  space  system  segment  may  be  considered  the  most 
important  since  the  other  two  require  the  space  segment  to 
be  operational  at  all  times.  Operationally,  the  space 
system  segment  consists  of  13  satellites  placed  in  orbits  at 
altitudes  of  10,398  nm  and  having  orbital  periods  of  about 
12  hours.   [Ref.  18] 
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The  original  plan  for  the  satellite  constellation 
consisted  of  24  satellites  in  three  planes  of  3  satellites 
each.  These  three  planes  were  to  be  placed  4S°  of  longitude 
apart  and  at  an  inclination  of  63  °  to  the  equator. 
[Ref.  19]  As  a  result  of  1980  budget  cutbacks  the  space 
system  segment  has  been  reduced  to  a  total  of  18  satellites, 
£R9f.  20]  Several  different  types  of  18  satellite 
constellations  have  been  proposed,  but  all  have  drawbacks. 
One  which  appears  the  most  promising  is  the  'Nonuniform  18* 
constellation  containing  6  planes,  each  of  which  contains  3 
satellites  spaced  120°  apart  (fig. 3.1.).  The  areas  of  user 
outages  (fig. 3. 2.)  are  a  disadvantage  of  this  configuration. 
The  user  outages  are  indicated  by  the  darker  areas  shown  on 
Figure  3.2.  These  outage  areas  will  appear  to  rotate  about 
the  globe  at  the  indicated  latitudes,  and  positioning  with 
the  GPS  will  be  intermittent,  or  of  decreased  accuracy  for 
users  at  these  latitudes.  The  final  constellation 
configurchtion  is  to  be  selected  early  in  1982.  [Ref.  2l] 

Sach  NAVSTAR  satellite  has  the  following 
characteristics:  weight  of  9  82  pounds,  length  of  17.5  feet 
from  tip  to  tip  of  the   5  square  solar  panels,   three  atomic 


44 


#   Proposed  Satellite 
O   jeleted  Satellite 


Fia.  3.1.   The  nonuniform  18  Satellite  Constellation 
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frequency  standards  (cloclcs)  ,  and  three  nickel-cadmium 
batteries  for  operations  in  darkness  and  during  peak  load 
periods.  Stabilization  of  each  NAVSTAR  is  maintained  by 
four  skewed  reaction  wheels  which  produce  three-axis 
stabilization.  NAVSTAR  station-keeping  and  momentum  damping 
capability  is  provided  by  a  hydrazine  propulsion  system 
within  each  satellite.  Life  expectancy  of  each  unit  is  five 
years,  and  expendables  are  expected  to  last  for  up  to  7 
years  with  normal  usage. 

Figure  3.3.  illustrates  the  parts  of  each  NAVSTAR.  A 
twelve  element  shaped-beam  helix  antenna  transmits  the 
navigation  message  on  two  L  band  frequencies  to  user  system 
segments.  Simultaneously  it  receives  satellite  status, 
clock  corrections,  ionospheric  data,  and  ephemeris  constants 
from  the  control  segment  on  two  S  band  frequencies. 
[Ref.  22]  The  two  coherent  L  band  frequencies  are  centered 
on  1227.6  MHz  (LI),  and  1575.4  MHz  (L2)  .  [ Ref .  23] 
C Ref •  24]  Both  frequencies  are  necessary  to  determine  the 
total  ionospheric  delay  correction  for  signal  propagation 
time  (Table  I  ) .  Errors  produced  by  the  ionospheric  delay 
in  signal  velocity  contribute  most  of  the  error  in  the 
pseudo-range.    By  utilizing   the  LI  and  L2   frequencies  and 
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Fig.    3.3.        NAVSTAR  Satellite 
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TABLE  I.   Ionospheric  Group  Delay  as  a  Function  of  Frequency 


^       -   R  ,  A  ,  B  ,  C 
GD 


C   ,2 


(neglecting  third  order  and  higher  terms  which  contribute  minimal 
error) 


where  A  =  ionospheric  delay  constant  of  conditions 

R  =  true  range 

C  =  speed  of  light 

f  =  carrier  frequency 

B  «>^  (averaged  earth  magnetic  field  strength  along 
the  path) 


GO 


=  group  delay  between  points  G  and  D 


Term  Healected; 
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Range  Error  at 
f  -  1  .5  GHz 

^  1  inch 
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the  following   algorithm,   the   error  is   decreased  but   not 
totally  eliminated.   [Ref,  25] 


T     ■  -  T    =   A  T 

GDL]   GDL2 


=  A 


Where: 


GDL 


^  GDL."   ^2 


"L. 


where  A  =  ionospheric  constant  and  "^  qq|_  ,  '^  qql  are  the 
respective  ionospheric  group  delays  at  frequencies  Li  and 
L2.  Both  frequencies  are  also  modulated  to  transmit  the  GPS 
time  codes  and  the  navigation  message. 

GPS  time  is  transmitted  using  two  pseudo-random  noise 
chip  codes — the  Coarse  Acquisition  (C/A)  and  the  Precision 
(P)  codes.  Modulation  on  the  L1  frequency  carries  both  the 
C/A  and  P  codes  in  phase  quadrature,  while  the  L2  frequency 
carries  only  the  P  code.  These  codes  identify  each  acquired 
satellite  by  matching  the  unique  pseudo-random  noise  code 
pattern  generated  by  each  satellite  with  similar  user 
generated  codes. 
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Both  codes  are  required  to  measure  the  navigation  signal 
propagation  time.  The  phase  shift  is  measured  and  used  to 
match  the  user  generated  signal  with  the  incoming  satellite 
signal,  [Bef.  26]  The  C/A  code  has  a  repetition  period  of 
one  millisecond  and  has  a  short  code  stream,  transmitting 
1.023  million  bits  per  second  (bps) .  This  code  is  easily 
acquired  and  matched  with  the  user  generated  pseudo-random 
noise  code  but  only  provides  a  coarse  time  signal  for 
computing  the  pseudo-range.  However,  the  P  code  has  a  long 
code  stream,  repeating  only  every  seven  days,  and  it  is 
transmitted  at  the  higher  rate  of  10.23  million  bps.  With 
this  higher  data  transmission  rate  and  long  repetition 
period,  a  more  accurate  time  signal  is  available  for 
pseudo-range  determination.  However,  the  long  period  and 
high  data  rate  makes  the  P  code  difficult  to  match  with  the 
user-generated  pseudo- random  noise  code  pattern  unless  the 
receiver  has  a  highly  accurate  time  standard,  and  the 
approximate  receiver  location  is  known.  In  lieu  of  this, 
the  usual  method  of  acquiring  the  P  code  is  by  matching  the 
C/A  code  of  the  desired  NA7STAR,  then  using  the  Handover 
Word  (HOW) ,  to  locate  the  correct  P  code  sequence  within  the 
7-day   pseudo-random      noise  code.  GPS   time      is   transmitted 
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with  the  HOW  every  six  seconds  and  provides  very  accurate 
GPS  time  to  the  user.  Both  of  these  codes  provide  the 
proper  time  reference  for  computation  of  the  pseudo-range. 
The  C/A  code  provides  a  lass  accurate  time  reference  than 
the  P  code.  The  selection  of  the  code  depends  only  on  the 
navigational  accuracy  desired. 

The  navigation  message  is  also  modulated  on  frequencies 
LI  and  L2  but  at  a  much  slower  rate  (only  50  bps)  than 
either  the  C/A  or  P  code.  On  both  frequencies,  the 
navigation  message  is  identical  and  consists  of  5  subsets  of 
6  seconds  duration  each.  Each  navigation  message  begins 
with  the  telemetry  message  (TLM)  and  the  HOW  message,  which 
are  used  to  transfer  from  the  C/A  to  ?  codes.  The  remainder 
of  the  entire  navigation  message  allocated  to  each 
subsection  is  transmitted  after  the  TLH  and  HOW  messages. 
Information  in  the  navigation  message  includes  approximate 
satellite  ephemerides,  status  of  other  NA7STAR  satellites, 
parameters  for  cloclc  corrections,  atmospheric  signal  delay 
corrections,  satellite  performance  status,  momentum  dump 
status  since  last  upload,  time  since  last  upload,  almanac 
and  identification  codes  for  all  satellites,  and  provisions 
for  inclusion  of  any  special   or  important  messages  from  the 
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control  system  segment  to  the  user.    [Hef.  27,  28J 
Navigation  is   thus  achievable  by  combining   the  information 
transmitted  by  the   space  system  segment  in  the  form  of  the 
C/A  and  P  codes  and  the  navigation  message. 

Accuracy  and  proper  function  of  GPS  are  governed  by  the 
second  segment — the  control  segment.  The  control  segment 
monitors  and  uploads  the  individual  NAVSTAR  satellites, 
checlcs  operational  characteristics,  determines  the  ephemeris 
and  almanac  data  and  corrections  reguired,  computes  the  time 
delay  and  clock  corrections  for  each  satellite  atomic 
standard,  determines  the  atmospheric  delay  corrections  for 
pseudo-range  computations  and  includes  special  messages  in 
the  body  of  the  navigational  message.  These  functions  are 
under  the  command  and  operational  control  of  the  Master 
Control  Station,  currently  located  at  Vandenberg  AFB, 
California.  Four  unmanned  remote  monitor  stations,  located 
in  Hawaii,  Guam,  Elmendorf  AFB,  Alaska,  and  Vandenberg  AFB, 
collect  the  pseudo-ranges  to  each  satellite,  the  change  in 
pseudo-range  of  each  satellite  signal,  local  meteorological 
data,  and  the  remote  monitor  station  atomic  standard 
parameters.  2ach  set  of  data  is  transmitted,  on  request,  to 
the  Master  Control  Station,   where   all  the  data  is  combined 
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and  the  satellite  upload  corrections  and  parameters  are 
computed.  These  corrections  and  parameters  are  then  relayed 
to  the  Opload  Station,  Vandenberg  AFB,  for  transmittal  to 
the  respective  NAVSTAR  satellites.  The  uploads  occur  when 
the  satellites  are  first  visible  to  the  Upload  Station  and 
when  the  User  Range  Error  exceeds  U  meters,  as  computed  by 
the   Master   Control  Station.     [Ref.  29,  30,  3lJ 

The  last  segment,  the  user  system  segment,  combines  the 
information  provided  to  the  space  system  segment  (by  the 
control  segment)  with  the  information  received  from  the 
space  system  segment  to  obtain  the  pseudo-ranges  from  four 
NAVSTAR  satellites.  These  ranges  are  computed  using  the 
following  algorithms: 

where  R^-  =  true  slant  range  to  satellite  i,    c  =  speed  of 

light  in  a  vacuum,   tt^  =  time  of  received  GPS  signal  by  user 

receiver  (assumed  to  be  simultaneous  from   all  satellites) , 

t^      =  time   of   transmission   at    satellite   i,     and 
i 

Aty^  =  atmospheric  propagation   delay   time  from  satellite  i 
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(fig.3.U.)«      Assuming  a    user   cloclc  error,    the    true    range   can 
be   expressed  in  terms  of   a  pseudo-range: 


R^      =    R.     +    cAt  „      ♦    c(At,-At^     ) 
P.  1  A.  ^^  R     ^  S. 

1  1  1 

where  Rp.  =  pseudo-range  including  clock  error,  ^""^r  ~  "tini© 
delay  in  the  user  receiver  of  the  received  signal,  and 
Ate  =  time  delay  in  satellite  transmission  from  satellite 
i. 


To  determine  the  slant   ranges: 


R=Ni{xc--x)     >   (yc  -y)  ^  ♦(Zc  -2)  ^ 

0.  :i.  b. 


R.      =^|(Xc    -X)^    ^{Yc    -y)^    >(Zc    -Z)^    +cAt,     ^c(At^-At       ) 

The  following  quantities  are  known:   x   ,  y   ,  z   obtained 

from   navigation  ephemeris  data;  A  t«    computed  using   the 

i 

ionospheric  delay  correction  obtained  from  observing 
frequencies  LI  and  L2  (tropospheric  delay  neglected)  ;  A  ts. 
obtained  from  the  navigation  message  as  the  satellite  clock 
error.  Thus,  there  are  four  unknowns:  x,  y,  z  (user 
position  coordinates)  and  ^  t  p^  (  user  clock  error  or  bias). 
To  solve  this  set  of  equations,  four  satellites  must  be 
observed  simultaneously.   [Ref.  32] 
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should  only  three  satellites  be  observable,  the  altitude 
may  be  assumed  and  the  solutions  can  be  computed  for  the 
other  two  location  variables  and  the  user  clock  error. 
Solutions  of  this  type  (altitude  hold)  significantly  degrade 
position  accuracy. 

The  above  description  is  over-simplified  since  the 
relativistic  effects  of  NAVSTAR  satellite  velocity,  as 
observed  by  the  user,  are  neglected  in  terms  of  the 
satellite  time  delay  correction  and  the  atmospheric 
propagation  delay  correction.  The  actual  user  velocity,  in 
terms  of  velocity  in  the  x,  y,  and  z  directions,  is 
calculated  using  the  Doppler  shifts  of  the  two  frequencies 
LI  and  L2  and  the  estimated  satellite  position  taken  from 
the  ephemeris  data.  With  the  relativisric  effects  as  an 
integral  part  of  the  data  set,  solution  of  the  navigation 
position  problem  becomes  more  complex.  Phase  I  receiver 
manufacturers  have  included  them  in  the  solution  of  the 
problem  by  several  different  methods.  Methods  of  solution 
depend  on  the  usage  of  the  receiver  (whether  static,  high  or 
low  dynamic)  and  the  expected  cost/benefit  ratio  for  the 
receiver  (in  terms  of  the  total  cost,  weight  and  size). 
Current  state-of-the-art  design  technology  and  navigational 
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Fig.  3.4.  Pseudo-Range  Description 
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algorithms,  developed  since  Phase  I,  will  be  included  in  the 
Phase  II  receivers  and  will  significantly  improve  accuracy 
of  future  receivers. 

B.   MANPACK  REC2IV2R 

The  Phase  I  receiver  used  for  this  research  was  the 
Texas  Instruments  GPS  Manpack,  also  known  as  the  Manpack 
Vehicular  User  Equipment  (MVUE) .  MVUE  is  a  singla  channel 
low  dynamic  receiver,  usable  for  receiver  velocities  from  0 
to  25  meters  per  second.  [Ref.  33]  System  configuration  is 
illustrated  by  figure  3.5.  Composing  the  system  are  the 
aanpack  control  data  unit  (CDU) ,  antenna  and  vehicle  mount, 
and  vehicle  power  converter.  The  Manpack  is  the  receiver 
unit  (fig-  3.6.)  with  provisions  to  connect  an  external  data 
output,  through  the  MVUE  Instrumentation  System  (MIS) 
interface,  figure  3.7.  Facilities  for  an  internal 
preamplified  antenna  mount,  or  an  external  preamplified 
antenna,  and  an  attachment  for  either  an  external  power 
source  or  the  integral  battery  pack  are  also  shown  (figure 
3.7).  Together,  the  MVOE  and  the  SBP/9900  microprocessor 
data  processing  unit  have  storage  capacity  for  47,360 
sixteen  bit  words.  This  storage  capacity  allows 
interrogation  of  NA7STAR  satellites   sequentially,   with  the 
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Fig.  3.7.   Cutaway  of  MVUE  Receiver 
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receiver  computing  a  position  fix  every  24  seconds.  The 
relativistic  velocity  effects  on  the  received  signals  are 
determined  by  matching  the  incoming  frequencies  with  the 
expected  frequencies.  The  user  velocity  relative  to  the 
satellite  is  determined  from  the  observed  frequency  shifts. 
The  nominal  frequency  shift  for  each  satellite  transmission 
is  computed  using  the  simultaneously  transmitted  ephemeris 
data.  Total  weight  of  this  receiver  unit  is  33.5  pounds, 
and  it  has  a  volume  of  17,958  cubic  centimeters. 

The  second  element  of  the  avOE,  the  CDU  (fig.  3.8.),  is 
the  input/output  device  for  the  non-automated  Manpack  user. 
Inputs  are  entered  by  keying  the  pressure  pads  set  in  six 
rows  of  four  pads  per  row.  To  obtain  a  navigational 
position,  the  approximate  user  latitude,  longitude,  and 
Greenwich  Mean  Time  (GMT)  are  keyed.  After  the  warmup 
period  of  3  to  5  minutes,  the  satellite  signals  are 
received,  and  the  CDU  will  display  the  actual  position. 
Additional  input  may  include  NA7STAR  satellite  selection, 
user  altitude,  and  other  position  datums.  Output  options 
include  position,  GMT  time,  and  system  error  messages. 
Other  inputs  and  outputs  are  listed  in  the  MVUE  manual. 
£Hef.  3U] 
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Fig.    3.3.        iMVUE   CDU   Keyboard/Display   Layout 
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Automated  users  may  elect  output  either  to  a  radio 
interface  or  through  the  MIS  interface.  For  this  research, 
a  special  interface  unit,  designed  and  constructed  at  the 
Office  of  Testing  and  Evaluation  of  Systems,  NOAA,  was  used 
to  link  the  MIS  interface  board  within  the  receiver  with  a 
Hewlett  Packard  9825T  computer.  With  this  interface,  user 
latitude,  longitude  and  GPS  time  were  logged  onto  magnetic 
tape  for  future  processing.  Figure  3.9.  and  Table  II  are 
the  engineering  documents  used  in  construction  of  the  MVOS 
interface  unit. 

Other  MVUE  elements  are  the  antenna  and  the  antenna 
external  preamplifier  vehicle  mount.  Located  within  the 
antenna  are  two  antenna  subassemblies  associated  with 
frequencies  LI  and  L2.  Both  of  the  subassemblies  are 
connected  to  the  external  preamplifier  with  two  lengths  of 
coaxial  cable,  enclosed  in  a  flexible  conduit.  The  entire 
unit,  antenna  and  preamplifier,  is  connected  to  the  Manpack 
receiver  with  a  power  cable  and  a  coaxial  antenna  cable. 
Cable  length  is  determined  by  the  location  of  the  external 
antenna  but  is  limited  by  the  minimum  strength  required  by 
the  Manpack  receiver.  The  combined  antenna  and  preamplifier 
unit  is  36  inches  long,  and  it  weighs  1.9  pounds. 


64 


u 


<u 

> 

<u 
u 

(U 

a: 


(T3 

S- 


3 

u 


65 


TABLE 

II.  MANPACK 

Interface 

Wiring  Code 

Wire  Color 

OUT 

IM 

Wire  Color 

Code 

BIT 
1 

1 

26 

BIT 
1 

Code 

90 

0 

91 

2 

2 

27 

2 

1 

92 

3 

3 

28 

3 

2 

93 

4 

4 

29 

4 

3 

94 

5 

5 

30 

5 

4 

95 

6 

6 

31 

6 

5 

96 

7 

7 

32 

7 

6 

97 

8 

8 

33 

8 

7 

934 

9 

9 

34 

9 

912 

935 

10 

10 

35 

10 

913 

936 

11 

11 

36 

n 

914 

937 

12 

12 

37 

12 

915 

945 

13 

13 

38 

13 

923 

946 

14 

14 

39 

14 

924 

947 

15 

15 

40 

15 

925 

948 

16 

16 

41 

16 

926 

98 

PCTL 

17 

42 

PFLG 

8 

927 

CTLO 

18 

43 

PSTS 

908 

928 

CTLl 

19 

44 

STIO 

915 

901 

1  BIT  1/0 

20 

45 

STIl 

917 

902 

PRESET 

21 
22 
23 
24 

46 
47 
48 
49 

EIR 

918 

GMD 

25 

50 

GND 

6b 


The  power  converter  must  supply  24  volts  DC  to  the 
Manpack  receiver.  During  this  research,  a  power  supply 
which  converted  120  volts  AC  to  24  volts  DC  was  used.  This 
power  supply  was  considered  more  stable  than  the  power 
converter  provided  with  the  MVOE,  and  the  voltage  output  and 
power  were  controlled  and  monitorad  during  the  testing 
period. 

C.   DIFFERENTiaL  GPS 

As  mentioned  earlier  in  this  chapter,  atmospheric 
propagation  delays  are  the  major  source  of  error  in  the 
navigational  equations  because  ionospheric  propagation  is 
poorly  understood.  The  tropo spheric  propagation  delays  are 
neglected  in  the  navigation  equations,  and  the  resultant 
error  is  dependent  on  local  temperature,  pressure,  humidity, 
and  NA7STAR  elevation  and  azimuth.  Tropospheric  corrections 
are  usually  neglected  if  the  observed  satellite  is  more  than 
5°  above  the  local  horizon.   [ Ref •  35] 

Additional  error  sources  are  receiver  and  satellite 
clock  perturbations,  lengths  of  signal  paths  from  different 
satellites,  multipathing  delays,  satellite  ephemeris  errors, 
and  receiver  noise.  Combined,  these  factors  lead  to  errors 
of  3.6  to  6.3  meters,  as  table  III  illustrates.    [Ref.  36] 
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TABLE  1 1 1. RANGE  ERROR  BUDGET 


Uncorrected  Error  Source 


User  Equivalent  Range  Error, 

(1°-)       ' 


feet 


meters 


SV  Clock  Errors 


Ephemeris  Errors 

5.0 

1.5 

Atmospheric  Delays 

8.0  -  17.0 

2.4  -  5.2 

Group  Delay  (SV  Equipment) 

3.3 

1.0 

Multipath 

4.0  -  9.0 

1.2  -  2.7 

Receiver  Noise  and  Resolution 
Vehicle  Dynamics 


5.0 


1.5 


Root  Sum  Square 


11.8  -  20.7 


3.6   -  6.3 


*Two  hours  after  update 
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One  method  used  to  reduce  the  magnitude  of  this  error  is 
that  cf  differential,  or  translocated,  positioning.  To  use 
this  method,  one  receiver  is  positioned  over  a  known 
geographic  position.  Pseudo-ranges  to  each  satellite  are 
determined,  and  a  differential  correction  is  obtained  for 
each  pseudo-range.  These  corrections  are  nhen  applied  to 
the  satellite  pseudo-range  observed  at  a  second  mobile 
receiver.  With  the  differential  corrections  applied,  the 
net  error  of  the  mobile  GPS  receiver  is  significantly 
reduced.   [Ref.  37] 
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IV.   TSST  PROCSDORES 

A.   POSITIONING  NETWORK 

One  of  the  considerations  in  planning  this  test  was  to 
examine  the  various  methods  of  positioning  the  shore-based 
equipment.  It  was  desirable  to  use  closely  spaced  (1-2  km) 
locations  within  the  immediate  surroundings  of  Monterey 
harbor  for  reasons  of  ready  access  by  personnel  during 
installation  of  equipment  and  for  optimum  accuracy  of  the 
positioning  systems  chosen  for  control.  As  mentioned 
previously,  the  positioning  systems  chosen  for  the 
comparison  (the  three  theodolites  and  the  ARTEMIS)  give  best 
results  at  relatively  short  ranges. 

Two  predetermined  station  marks  were  selected  for  the 
test — 1)  ass  HON  1978,  a  bronze  disk  set  by  the  Array  Corps 
of  Engineers,  located  on  a  sand  dune  on  the  Naval 
Postgraduate  School  property  4  00  meters  from  the  ocean  at  an 
elevation  of  ao  meters  above  sea  level  and  2)  MOSSEL  1932,  a 
bronze  disk  set  by  the  Coast  and  Geodetic  Survey  in  a  rock 
outcrop  7  meters  above  sea  level  at  the  shoreline  on  the 
property  of  the  Hopkins  Marine  Laboratory  of  Stanford 
University.   The  third  station  mark  was  set  in  the  center  of 
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the  Monterey  Coast  Guard  Pier  for  purposes  of  the  test. 
This  was  a  DMA  bronze  disk  stamped  SE0CEI7EH  STA  31370. 
This  location  was  chosen  for  access  to  power  facilities 
available  at  the  pier  and  because  of  central  location  to  the 
test  area. 

The  position  of  the  Geoceiver  station  was  established 
using  a  Magnavox  Geoceiver  (Model  AN/PRS  14)  obtained  from 
the  Department  of  Satellite  Geophysics,  Satellite  Tracking 
Branch,  DMA/Hydrographic  Topographic  Center.  The  equipment 
was  operated  in  accordance  with  DMA  requirements  to  obtain  a 
documented  first  order  station.  The  data  were  submitted  to, 
and  processed  by,  DMA/HTC  in  Brookmont,  Maryland,  using 
precise  ephemeris  information  to  provide  the  WGS-72 
position.  Michael  211ett  of  DMA  provided  the  technical 
expertise  for  acquisition  of  the  data  and  maintenance  of  the 
gecceiver  during  the  test  period. 

A  third  order.  Class  I,  closed  traverse  was  conducted 
following  the  standards  of  the  National  Geodetic  Survey 
(NGS) ,  using  a  Wild  T-2  theodolite  to  measure  angles  and  a 
Motorola  MRA  V  Tellurometer  to  measure  distances.  The 
THANSIT  position  computed  by  DMA/HTC  was  used  for  the 
initial  position  of   the  survey,   and  a   Polaris  observation 
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using     the  theodolite  provided     the      initial   azimuth       (fig. 

4.1.).        A    total      closure  error  in   distance      of   0.061    meters 

and   angle   of      5.U"    was  achieved.        The      positions   determined 

from   the      traverse  and  used   for     the  control    networic      are   as 

follows: 

OSE    aON  36^  36«08  .a24"N            121°  52»  40 .  301"W 

MUSSEL  36°  37»21  .832"N            1 2  f  54  •  16  .  048"W 

GEOCEIVEfi    STA    31370  36°  36«  36  .246"N            121°  53*  29  .693"W 

B.       POSITIONING    EQUIPMENT 

As  previously  described,  three  positioning  systems  were 
used  during  the  test — the  three  theodolites^  the  GPS,  and 
the  ARTEMIS.  In  addition,  the  geoceiver  was  set  to  acquire 
data  during  the  shipboard  portion  of  the  test.  During  the 
time  of  day  that  the  GPS  satellites  were  available  (Table 
IV)  only  one  or  two  observations  per  day  of  the  same 
THANSIT  satellite  were  possible.  The  geoceiver  data  were 
recorded  but   were  not  used   for   this   study. 
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Fig  .  4.1.   Test  Area  with  Traverse  Sketch 
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1 .   Theodolites 

This  section  describes  the  placement  and  operation 
of  the  equipment  during  the  test.  One  theodolite  was  set  on 
a  tripod  at  each  of  the  three  stations,  and  each  was 
occupied  by  an  observer  and  a  recorder.  Communication 
between  stations  was  maintained  via  portable  radios. 

Because  tests  could  be  conducted  only  between  10  PM 
and  2  AM  local  mean  time  (LMT)  due  to  GPS  availability 
times,  it  was  necessary  to  plan  a  systems  of  lights  for  the 
stations.  The  three  theodolite  sites  were  made  visible  by 
placing  battery  powered  lights  directly  over  the  station 
marks.  In  addition,  battery  operated  light  packs  for  the 
theodolites  were  used.  All  shorebased  lights  and  light 
pacJcs  for  the  test  were  provided  by  the  Pacific  Marine 
Center,  NOAA. 

The  ship  system  was  designed  with  an  omnidirectional 
pair  of  lights  anchored  around  the  support  pipe  for  the  GPS 
antenna.  Also,  the  equipment  shelter  on  the  boat  deck  of 
the  ship  was  lighted  by  spotlights  so  that  the  GPS  antenna 
was  easily  seen  by  the  shorebased  theodolite  operators. 

During  the  test,  each  theodolita  observer  initialled 
on  one  of  the  other  stations  and  recorded  the  reading  on  a 
log  sheet.     The  theodolita  operators  then   intersected  the 
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lights  on  the  ship  GPS  antenna  using  the  theodolite 
telescope  and  tracked  the  ship  until  a  'aark«  was  called 
over  the  radio  from  the  personnel  aboard  the  ship.  The 
angle  measured  at  this  time  was  recorded  on  the  log  sheet  as 
a  fix.  The  time  interval  between  fixes  was  approximately 
one  minute.  2ach  series  of  approximately  twenty  fixes 
constituted  a  set.  At  the  end  of  each  set,  a  final  reading 
was  made  of  the  angle  to  the  reference  station.  The  fixes 
were  numbered  consecutively  throughout  the  test,  and  the 
Greenwich  Mean  Time  (GMT)  of  each  fix  was  recorded  on  the 
ship  log  sheets.  The  length  of  the  sets  depended  on  the 
ship  position  and  heading. 

To  correct  the  observed  theodolite  angles,  the  mean 
of  the  final  pointing  and  the  initial  pointing  was 
determined  for  each  set  and  for  each  station.  This  value 
was  assumed  to  be  a  constant  throughout  the  set  and  was 
removed  from  each  angle.  A  collimation  error  for  each 
theodolite  was  also  removed  from  each  theodolite 
observation.  The  same  observers  and  theodolites  were  used 
nightly  at  all  stations  except  MOSSEL,  where  two  different 
individuals  observed  over  the  test  week.  This  continuity  of 
personnel  and  equipment  was  designed  to  reduce  operator  and 
instrument  error  to  a  minimum. 
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2.   GPS  and  ARTS  MIS 

To  colocate  the  different  pieces  of  equipment  at  the 
GEOCEIVER  STA  31370,  a  wooden  tower  was  built  as  shown  in 
figure  4.2.  This  tower  was  set  over  the  station  for  the 
period  of  the  test.  An  NPS  equipment  vehicle  was  parked  as 
close  as  possible  to  the  tower.  Electrical  power  was 
obtained  from  the  Monterey  Coast  Guard  facility. 

Prior  to  the  ship  test,  both  GPS  antennas  were 
placed  on  the  tower  and  centered  over  the  station  mark.  ?or 
the  first  two  nights,  9  and  10  Hay,  equipment  problems 
allowed  the  use  of  only  one  system  each  night  and  data  were 
acquired  only  to  determine  the  stability  of  each  system. 
During  the  next  night,  11  May,  position  outputs  for  the  two 
simultaneously  operating  GPS  receivers  were  recorded.  This 
static  test  gave  a  determination  of  any  difference  between 
the  two  receivers  under  almost  identical  conditions. 

For  the  ship  test,  the  GPS  antenna  was  centered  over 
the  mark  and  was  anchored  at  the  top  of  the  tower;  the 
geoceiver  antenna  was  centered  over  the  mark  and  anchored  at 
the  middle  level;  and  the  theodolite  on  the  tripod  was  set 
over  the  mark  at  ground  level.  This  arrangement  provided 
colocation  of   all  three  systems.     It  was  assumed   that  no 
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error  was  introduced  by  the  presence  of  the  wooden 
structure.  The  shorebased  ARTEMIS  antenna  unit  was  offset 
from  the  theodolite  to  eliminate  interference  in  the  signal 
due  to  the  tower.  An  offset  position  was  computed  relative 
to  the  station  mark  and  was  used  in  postprocessing  of  the 
ARTEMIS  data. 

Except  for  the  antennas  for  the  GPS,  the  ARTEMIS, 
and  the  geoceiver  systems,  all  other  equipment  was  installed 
in  the  equipment  vehicle,  where  an  operator  was  stationed 
throughout  the  test  to  monitor  the  data  quality. 

The  GPS  Manpack  is  a  manually  operated  system.  As 
mentioned  previously,  it  was  decided  that  manual  operation 
of  GPS  was  detrimental  to  the  goal  of  the  test,  and  that  an 
interface  that  would  provide  an  automated  printout  of 
positions  from  the  Manpack  as  well  as  a  printout  of  GPS  time 
would  be  desirable.  Knute  Berstis  and  Gary  Hhitsell  of  the 
Office  of  Testing  and  Evaluation  of  Systems  (OTES) ,  NOAA, 
assembled  two  interfaces  that  provided  data  storage  on 
cassette  tapes  as  well  as  the  desired  printout.  Each  system 
operated  in  conjunction  with  a  HP-9825  computer  and  a 
HP-78a5  printer/plotter.  Software  was  provided  by  OTES, 
NOAA.    This  arrangement  permitted  storage  of  the  satellite 
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messages  for  further  processing  and  provided  an  automated 
printout  of  the  ARTEMIS  data.  A  GPS  ManpacJc  had  been  sent 
to  OTES  for  testing  during  the  construction  of  the 
interface.  Technical  expertise  and  operation  of  the 
interface  during  the  test  were  provided  by  OTES. 

The  data  from  the  GPS  receiver  were  printed  at  an 
interval  of  24  seconds,  and  consecutive  data  were 
continually  recorded  on  tape  and  printed  during  the  test 
periods.  In  addition,  a  manual  record  of  the  GPS  position 
from  the  CDO  was  logged  as  the  • marks*  were  called  via  the 
radios. 

The  ship  positioning  systems  were  installed  as  shown 
in  figures  4.3.  and  4,4.  An  equipment  shelter  was  secured 
to  the  boat  deck  of  the  R/V  ACANIA.  A  length  of  six-inch 
diameter  polyvinyl  chloride  (PVC)  pipe  was  attached  to  the 
aft  wall  of  the  building  on  the  centerline  of  the  ship.  The 
GPS  antenna  and  the  system  of  lights  previously  mentioned 
were  set  on  top  of  the  pipe. 

The  ARTEMIS  mobile  antenna  unit  was  set  on  a 
platform  located  on  the  centerline  of  the  ship  and  at  the 
top  of  the  building  high  enough  to  not  interfere  with  the 
signal  from   the  ship   radar  beam.    For   the  period   of  the 
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Fig.  4.4.   Equipment  Installation 
on  the  R/V  ACANIA 
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test,  the  weather  conditions  allowed  ship  operations  without 
use  of  the  radar,  thus  there  was  no  possibility  of  damage  to 
the  ABTEMIS  unit. 

All  other  equipment — the  ARTEMIS  receiver  unit,  the 
GPS  receiver,  control  data  unit,  power  filter,  and 
interface,  the  HP-9825,  the  HP-7845  plotter/  printer,  and 
the  real  time  clock — was  placed  inside  the  building.  An 
operator  monitored  the  equipment  during  the  ship  tests  and 
another  individual  radioed  the  'mark*,  or  fix,  to  the 
theodolite  operators  on  shore.  The  equipment  used  at  the 
shore  location  was  installed  in  a  vehicle  and  consisted  of 
the  same  components  as  on  the  ship.  Operation  of  equipment 
during  the  ship  tests  was  monitored  by  an  individual  located 
in  the  vehicle. 

C.   LOG  OF  EVENTS 

The  test  was  divided  into  two  portions — 1)  a  shorebased 
test  which  consisted  of  two  nights  of  stability  testing,  and 
one  night  of  receiver  comparison  tests,  and  2)  the  ship  test 
consisting  of  three  nights  of  shipboard  operation.  The  R/V 
ACANIA  from  the  Naval  Postgraduate  School  was  the  platform 
used  for  the  shipboard  tests. 
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1.  stability  Xs^ 

Operations  began  on  9  May,  1981,  GMT.  One  of  the 
GPS  systems  was  placed  over  the  station  mark,  and  continuous 
data  was  recorded  during  the  satellite  availability  periods. 
In  this  manner,  system  1  was  tested  on  9  May,  and  system  2 
was  tested  on  10  May.  A  position  was  defined  as  a  quality 
observation  if  four  or  more  satellites  were  used  by  the 
receiver  to  obtain  the  output. 

2.  Static  Test 

Both  GPS  antennas  were  centered  over  GEOCEIVER  STA 
31370,  and  data  from  simultaneous  operation  of  both 
receivers  were  recorded  at  the  equipment  vehicle.  One 
antenna  was  placed  on  top  of  the  wooden  structure  and  the 
other  immediately  below  it  on  the  middle  level.  The 
vertical  displacement  was  assumed  to  cause  negligible  error 
in  the  horizontal  positions.  Data  were  acquired  on  1 1  May 
from  both  receivers  with  a  quality  observation  defined  as  a 
position  determined  from  four  or  more  satellites. 

3.  shipbgflgg  -gst 

?or  the  entire  period  of  shipboard  testing,  the  wind 
speeds  were  less  than  5  knots  and  the  seas  were  calm  with  1- 
2  feet  or  less  of  swell.  Visibility  was  excellent  throughout 
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the  test  week.  The  ship  track  for  each  set  is  shown  in 
figure  4.5.    A  log  of  subsequent  operations  follows: 

12  May The  GPS  Manpack  receiver  was  installed  on  the  R/V 

ACAHIA  and  SPS  and  theodolite  operations  were  conducted. 
The  ship  was  allowed  to  drift  during  sets  1-3  (fix  1-69) 
with  23  fixes  per  set,  in  a  east-northeasterly  direction. 

13  May The  ship  cruised  at  U  knots  on  track  lines  from 

south  to  north  for  sets  4  and  6,  and  lines  from  north  to 
south  for  sets  5  and  7.  During  fixes  139  through  157  the 
radio  at  the  Coast  Guard  pier  was  disabled,  and  no  data  were 
acquired  by  either  the  theodolite  observer  or  the  operator 
in  the  equipment  vehicle. 

14  May The  shipboard  GPS   system  was  inoperable.     The 

trouble  was  traced  to  a  faulty  antenna  cable,  requiring 
replacement.   No  data  were  acquired, 

15  May The   ARTEMIS,   the  GPS   and  the  theodolite  network 

were  all  used  for  this  test  period.  Data  from  sets  9 
through  13  were  acquired  (fixes  165-262)  . 

16  May The  control  data  unit  on  the  shore  based  GPS  system 

failed  before  the  start  of  ship  operation  and  could  not  be 
repaired  without  factory  maintenance. 
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k  total  of  2  52  fixes  v#ere  recorded  by  the  theodolite 
network  and  all  GPS  positions  corresponding  to  the 
theodolite  fixes  were  manually  logged.  The  ARTEMIS  azimuth 
and  distance  information  and  consecutive  GPS  positions  on 
the  ship  and  in  the  vehicle  were  stored  on  cassette  tapes 
and  printed  by  the  HP-7845  printer/plotter. 
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7.  DATA  PROCESSING  M2TH0DS 
Processing  of  the  thesis  data  was  divided  into  four 
interrelated,  but  separate,  tasks.  These  tasks  were  the 
basic  geodetic  survey,  theodolite  and  ARTEMIS  ship  position 
determination,  differential  GPS  ship  position  derernination, 
and  the  statistical  analysis  of  the  ship  position 
ccoparisons.   Each  task  will  be  discussed  separately. 

A.   BASIC  SORVEI  COMPOTATIONS 

As  described  earlier  in  the  Test  Procedures,  a  basic 
geodetic  survey  of  the  test  area  was  conducted  to  place  the 
theodolite  control  stations  on  the  World  Geodetic  System  of 
1972,  or  WGS-72.  To  begin  the  survey,  a  aagnavox  AN/PRR  14 
Geoceiver  was  placed  over  the  center  survey  disk,  GEOCEIVER 
STA  31370  1981,  and  the  digital  Doppler  data,  from  TRANSIT 
satellite  68,  were  recorded  onto  paper  tape.  Thirty-three 
passes,  over  a  period  of  two  and  one  half  weeks,  were 
recorded,  but  analysis  of  the  equipment  maintenance 
diagnostics  indicated  a  possible  error  in  the  records. 
These  data  were  therefore  discarded. 

A  second  set  of  forty  satellite  passes  of  TRANSIT 
satellite  63  were  recorded  over  a  ten  day  period  and  were 
determined  to  be  of  acceptable   quality  according  to   the 
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standards  of  the  DMA  Department  of  Geodssy.  This  second  set 
of  data  was  then  sabaitted  to  DMA,  Department  of  Geodesy, 
Computation  Section,  for  analysis.  Also  submitted  were 
meteorological  data  taken  during  the  passes  and  additional 
information  required  for  documentation  of  a  first  order 
TRANSIT  position.  Satellite  68  Doppler  data  for  the  survey 
was  then  processed  in  conjunction  with  the  tracked  satellite 
position  by  DMA  computer  facilities.  A  final  geodetic 
position  on  the  WGS-72  ellipsoid  was  obtained  for  GEOCEIVER 
STA  31370.  Coordinates  of  the  station  are  listed  in  the 
Test  Procedures  section.  The  computed  height  relative  to 
the  ellipsoidal  surface  is  -33.39  meters.  Computation  of 
the  geographic  position  utilized  a  minimum  of  38  individual 
TRANSIT  satellite  passes  recorded  for  GEOCEIVER  STA  3  1370. 
This  provided  for  an  accuracy  of  one  part  in  10^  (first 
order  accuracy  standard).  [Ref,  38] 

To  orient  the  positioning  network,  an  initial  azimuth 
from  GEOCEIVER  STA  3  1370  to  station  MUSSEL  1932  was  computed 
from  two  sets  of  observations  of  POLARIS.  All  observations 
and  computations  of  the  star  sighting  utilized  procedures 
outlined  in  the  National  Ocean  Survey  Photogrammetric 
Instructions  Number  4,  Revision  1,  U/24/73,  the  Hydrographic 
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Manaal,   and  the  Manual  of  Geodetic  Triangulation,   SP-2U7. 
[Ref.  39,  403 

Having  the  initial  azimuth,  distances  were  measured 
between  all  three  control  stations  using  a  Telluroneter  MRA 
5  microwave  distance  measuring  unit.  Twenty  observations  of 
distance  were  made  for  each  of  the  three  lines  measured: 
GEOCEIVEH  STA  31370  to  MOSSEL ,  GE0CEI7EB  STA  31370  to  OSE 
MONr  and  MUSSEL  to  USE  MON.  The  average  of  the  distance  for 
each  line  was  then  corrected  for  instrument  error  and 
meteorological  conditions  following  the  instructions 
contained  in  the  Operators  Manual  and  the  Tellurometer 
Manual.  [Ref.  41]  Resulting  slope  distances  were  then 
reduced  to  geodetic  distances  using  the  respective 
corrections  for  slope  to  horizontal  distance,  reduction  of 
elevation  to  sea  level,  and  the  chord  to  arc  distance. 
These  corrections  were  computed  using  the  WGS  1972 
ellipsoidal  parameters  and  the  methods  outlined  in  the 
Hydrographic  Manual,  the  Manual  of  Geodetic  Triangulation, 
and  Introduction  to  Surveying.  [Ref.  U2]  The  resulting 
geodetic  distances  meet  third-order,  class  I  specifications 
for  traverse  accuracy,  1:20,000.  In  addition,  angles  to  the 
other  stations  were  observed  at  each  control  station,  giving 
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a  total  of  three  sets  of  angles.  Each  set  was  comprised  of 
four  direct  and  four  reverse  plate  observations  with  a  Wild 
T-2  theodolite.  Each  pair  of  direct  and  reverse 
observations  was  corrected  for  initial  pointing  errors. 
Observations  and  computations  were  made  using  the  methods 
detailed  in  the  Hydrographic  Manual,  the  Manual  of  Geodetic 
Triangulation,  and  the  Manual  of  Second  and  Third-Order 
Triangulation  and  Traverse.  [Bef.  43]  The  resulting 
closure  meets  the  specifications  for  traverse  accuracy  of 
three  seconds  per  traverse  station  or  a  closure  of  6 
seconds. 

Hith  the  geodetic  position  of  GEOCEIVEH  STA  31370  known, 
as  well  as  the  azimuth  from  GEOCEIVER  STA  31370  to  HOSSEL 
and  the  distance  from  GEOCEIVER  STA  31370  to  MOSSEL,  the 
geographic  position  of  MUSSEL  was  computed  on  the  iGS  1972 
ellipsoid  using  the  direct  geodetic  computation.  This 
computation  used  the  Rainsford's  method  modified  by  T. 
Vincenty  and  including  the  Helmert's  elliptical  terms. 
[Hef.  44] 

Next,  using  the  azimuth  from  MUSSEL  to  OSE  MON,  the 
geodetic  distance,  and  the  results  of  the  direct  geodetic 
computations,    the   geographic  position   of   USE   MON   was 
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determined.  Also,  the  ARTEMIS  geographic  position  was 
calculated  as  an  offset  using  an  azimuth  and  distance  from 
GEOCSIVER  STA  31370.  With  these  computations,  the  geodetic 
survey  vas  completed. 

Positioning  of  the  ship,  R/V  ACANIA,  required  one 
azimuth,  or  line-  of-position,  from  each  of  the  control 
stations  MUSSEL,  GEOCEIVER  STA  31370,  and  USE  MON.  These 
lines  of  position  were  used  to  triangulate  the  vessel 
position.  Angles  from  each  of  the  three  theodolite 
locations  were  observed  as  described  in  the  Test  Procedures. 
Each  of  the  three  angular  measurements  were  then  corrected 
for  mean  pointing  errors  and  the  individual  theodolite 
collimation  errors.  Each  corrected  angle  was  then  adjusted, 
by  control  station,  to  obtain  an  azimuth  from  south,  as 
required  to  triangulate  the  vessel  position.  The  three 
resulting  azimuths  and  the  corresponding  control  station 
latitude  and  longitude  were  then  used  in  a  least  squares 
three  line  of  position  intersection  computation.  Results  of 
the  computation  were  a  vessel  position  in  latitude  and 
longitude,  error  ellipse  semi-major  radius,  error  ellipse 
semi-minor  radius,  azimuth  from  south  of  semi-major  radius, 
and  circular  standard  error.  [Ref.  *5 ]  Figure  5.1 
illustrates  these  error  ellipse  parameters. 
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THE  ERROR  ELLIPSE 


The  mean  displacement  of  a  point  to  be  expected  in  the  direction  of 
the  two  coordinate  axes,  X  and  Y,  due  to  measuring  errors,  can  be  ex- 
pressed by  the  quantities  t  and  cr  . 

X      y 

When  the  north-south  axes  are  rotated,  two  axes  (X  and  Y)  perpendicular 
to  one  another  are  obtained,  corresponding  to  the  maximum  and  minimum 
values  of  ^   and  cr  ,  The  maximum  and  minimum  values  of  a-,  and  the 
clockwise  rotation  aiigle  (<=><)  of  the  major  axes  of  the  error  ellipse  are 
included  as  part  of  the  output  of  the  Geodetic  Package  and  can  be  used 
as  input  if  desired. 

The  circular  standard  error  (CSE)  is  also  included  in  the  output  with  the 
error  ellipse  parameters.  The  value  of  CES  is  the  radius  of  a  circle 
whose  area  is  equal  to  the  area  of  the  error  ellipse.  That  is: 


CSF  = 


max 


min 


EXAMPLE:   Error  Ellipse 
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Fig.    5.1 
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For  this  study,  262  vessel  positions,  or  fixes,  were 
observed  and  recorded.  On  examination  daring  processing,  51 
fixes  were  found  to  contain  errors  and  data  omissions  which 
rendered  the  vessel  positioning  computation  invalid  or 
indeterminate.  The  remaining  211  fixes  comprised  the  data 
set  which  was  then  broken  down  into  smaller  sets  for 
analysis  and  comparison.  These  sets  were  numbered  1  through 
13  and  categorized  according  to  vessel  trade.  All  sets  were 
then  processed  using  several  different  methods  to  determine 
the  most  accurate  approach. 

B.   TBEODOLITE  AND  ARTEMIS  SHIP  POSITION  DETERMINATION 

Before  processing  the  GPS  Manpack  data,  the  ARTEMIS  data 
were  processed  to  ensure  the  vessel  positions,  as  obtained 
by  theodolite,  were  correlated  with  the  correct  automated 
and  hand  recorded  GPS  data.  ARTEMIS  data  consisted  of  an 
angle,  MOSSEL-ARTEMIS  OFFSET- vessel,  and  distance  from  the 
ARTEMIS  OFFSET  to  the  vessel.  The  angular  measurement  was 
first  converted  to  an  azimuth  from  south  using  the  azimuth 
of  ARTEMIS  OFFSET  to  MUSSEL  as  the  correction  and  adding  the 
measured  angle.  Vessel  position  based  on  ARTEMIS  was  then 
computed  using  the  direct  geodetic  program  (also  used  for 
the  basic  survey) .    Comparisons  of  the  resulting  geographic 
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positions  and  the  theodolite  derived  geographic  positions, 
were  lade  for  corresponding  fix  numbers  using  the  inverse 
geodetic  computation  on  the  WGS-72  allipsoid.  The  inverse 
computation  utilizes  two  geographic  positions  as  inputs  for 
calculating  the  forward  and  back  azimuths  and  the  distance 
between  the  positions.  [Hef,  46]  For  comparison,  if  the 
two  positions  differed  by  more  than  20  meters,  the  ARTEMIS 
positions  on  either  side  of  the  suspect  ARTEMIS  fix  were 
then  used  to  compute  the  inverse  computation  to  the 
theodolite  ship  position.  In  seven  fixes,  the  fix  position 
values  were  changed  in  the  May  15  GPS  data  to  give  more 
accurate  results. 

C.   DIFFERENTIAL  GPS  SHIP  POSITION  DETERMINATION 

Hith  the  check  and  adjustment  of  the  GPS  data  and  the 
theodolite  determined  vessel  position  complete,  the  first 
GPS  data  processing  was  the  comparison  of  the  uncorrected, 
non-differential  ship  GPS  Manpack  receiver  data  to  the 
theodolite  computed  ship  position.  Using  the  inverse 
geodetic  computation,  the  forward  azimuth  and  the  distance 
between  the  theodolite  computed  ship  position  and  the 
ncn-dif ferential  GPS  ship  position  were  completed.  These 
results  were  then   used  in  the  analysis  of  the  differential 
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processing  to  examine  the  stability  of  the  shipboard  Hanpack 
receiver  compared  to  the  GPS  Manpacic  receiver  located  over 
GE0CEI7ER  STA  31370. 

There  were  two  methods  of  analysis  of  the  differential 
GPS.  The  first  of  these  was  the  determination  of  a  static 
correction  to  be  applied  to  all  the  GPS  data.  This  static 
correction  in  the  differential  GPS  computations  was  to 
account  for  a  suspected  bias  between  the  two  similar,  but 
not  identical,  Manpack  receivers.  A  mean  static  correction 
was  computed  for  the  most  stable  reception  period  during  the 
static  tests  on  11  May.  With  the  two  receivers  colocated 
over  the  GEOCEIVER  STA  31370,  a  mean  correction  of  -0.954 
seconds  of  latitude  and  -0.3  07  seconds  of  longitude  was 
found  between  GPS  Hanpack  System  1  and  System  2  during  times 
07:16:00  GMT  and  07:25:12  GMT.  This  static  correction  was 
then  added  to  the  latitude  and  longitude  of  GPS  Manpack 
system  1  (located  over  GEOCEIVER  STA  31370)  data. 

The  second  method  of  analysis  or  processing  in  the 
differential  mode  used  the  latitude  and  longitude  obtained 
directly  from  the  GPS  Manpack  receivers.  la  either 
processing  method,  the  computation  of  the  differential 
corrector  to  the  shipboard  GPS  Manpack  receiver  data  was  the 
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same.  An  explanation  of  the  computations  and  processing  of 
the  differential  correction  for  both  methods  follows. 

Computation  of  the  differential  corrector  for  each  GPS 
fix  position  was  accomplished  by  using  the  geographic 
position  obtained  with  the  stationary  Hanpack  receiver, 
located  over  GEOCEIVER  STA  31370,  and  the  geographic 
position  of  GEOCEIVER  STA  31370.  These  two  positions  were 
then  used  ia  the  inverse  geodetic  computation  to  obtain  a 
differential  correction  vector  (forward  azimuth,  and 
distance)  from  the  stationary  GPS  receiver  to  GEOCEIVER  STA 
31370,  In  figure  5.2,  the  differential  correction  vector  is 
labelled  "A".  Using  the  direct  geodetic  computation,  the 
differential  correction  vector  was  applied  to  the 
corresponding  shipboard  GPS  Hanpack  receiver  position,  and 
the  differential  GPS  position  was  obtained.  These  resulting 
positions  were  computed  on  a  one-to-one  basis  using  the  fix 
numbers  as  identifiers.  Vector  "3",  the  differential 
corrector  vector  which  is  identical  to  the  vector  "A"  in 
figure  5,2,  is  applied  to  the  shipboard  GPS  Manpack  receiver 
and  results  in  the  differential  GPS  ship  position. 

Comparison  of  the  differential  GPS  ship  position  with 
the  theodolite-computed  ship  position  was  accomplished  using 
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the  inverse  geodetic  computation  for  both  methods  of 
differential  GPS  processing.  The  resulting  forward  azimuth 
and  distances  for  both  processing  methods  were  used  in  the 
statistical  analysis. 

D.   STATISTICAL  DATA  ANALYSIS 

Statistical  analysis  was  done  on  data  separated  into 
sets  according  to  the  method  of  data  processing.  The  static 
data  and  stability  data  were  graphed  as  functions  of 
latitude  and  longitude  versus  time,  Oser  range  error, 
furnished  by  the  Space  Division,  GPS/OL-AQ,  Vandenburg,  AFB, 
was  plotted  as  a  function  of  time  for  each  day  of  the  test 
(see  Results).  Other  graphed  relationships  are  the 
geographic  positions  of  the  stationary  GPS  Manpack  receiver 
data  versus  the  GEOCEIVER  STA  31370  and  the  theodolite 
semi-aajor  radius  as  a  function  of  fix  number  (see  Results) . 

Other  statistical  data  analyses  were  represented  by 
histograms.  Plotted  on  these  histograms  are  the  relative 
frequency  on  the  left  margin,  variable  bar  increments  scaled 
on  the  bottom,  number  of  variable  observations  in  each  bar 
increment  on  the  top,  and  the  total  sample  size  on  the  top 
margin.  Within  the  histogram  the  frequency  bars  are  denoted 
by  asterisks,    the  mean  by  a   capital  M,   the   quartiles  by 
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periods  and  the  empirical  density      function   by   the   capital   F 
curve.      The  empirical   density    function   was  computed   using: 


1  N  X-Z 

F(Z)     =    Z    W      — J 

N(B(N))     i=1  B  (N) 


[B€£.  47]  where:  F(z)  =  the  empirical  density  function,  z 
=  the  frequency  of  observation,  x  =  individual  observation, 
N  =  the  total  number  of  observations,  H  =  weight  function,  B 
=  bandwidth, 

ff  =  1  -  z 
B(N)  =  range  of  all  observations  /  W 
The  statistical  values  listed  below  the  histogram  are 
defined  in  Table  7.  Computed  histograms  show  the  distances 
between:  1)  the  theodolite-computed  ship  position  and  the 
differential  GPS  ship  position,  both  with  and  without  the 
static  correction,  by  set;  2)  the  theodolite-computed  ship 
position  and  the  differential  GPS  ship  position,  without 
static  correction,  for  five  satellite  observation  periods  on 
each  of  the  three  ship  test  days;  3)  the  theodolite-computed 
ship  position  and  the  non-differential  GPS  ship  position  for 
the  five  satellite-observation  periods  on   each  of  the  three 
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TABLE  V.  STATISTICAL  DEFINITIONS 

The  definitions  of  the  basic  statistics  computed  by  the  program  are 
listed  below: . 

MEAN:   Average  of  the  sample 

MEDIAN:  Mid-Value  of  the  sample,  if  there  are  an  odd  number  of 
sample  points,  or  the  average  of  the  two  middle  values 
for  an  even  number  of  points 

TRIMEAN:   0.25  *  (Q^  +  2Q2  +  Q3),  where  Q^.  are  the  quartiles 

MIDMEAN:   The  average  of  all  sample  values  lying  between  the  upper 
and  lower  quartiles 

MIDRANGE:   Average  of  the  maximum  and  minimum 

GEOMETRIC  MEAN:   GM  =^>ix,'X^-.    .    .  x„ 

1  „2  n 

HARMONIC  MEAN:   HM  =  n  7^4] ^/^i 

VARIANCE:   Unbiased  estimators  for  variance  and  standard  deviation 
are  used  -  the  square  of  the  standard  deviation 

STANDARD 


DEVIATION:   S  =j^x^  -  n^ 

V  ' — 1 = — 

'    n  -  1 


COEFFICIENT  OF  VARIATION:   Standard  deviation  /  |Mean|  when  the 

mean  is  less  than  10.  **  -30,  the  coeffic- 
ient of  variation  is  set  to  zero. 

MEAN  DEVIATION:   The  average  of  the  sum  of  the  absolute  differences 

between  the  sample  values  and  the  median 

RANGE:   Maximum  -  minimum 

MIDSPREAD:   Q^  -  Q^ ,  also  called  the  interquartile  distance 

M3:   Third  central  moment.  Unbiased  estimator  is  used 

M4:   Fourth  central  moment.  Unbiased  estimator  is  used 

COEFFICIENT  OF  SKEWNESS:  MS  /  (Standard  Deviation)  **  3 

COEFFICIENT  OF  KURTOSIS:   (M4  /  (Standard  Deviation)  **  4)  -  3 

BETAl :   Biased  estimate  of  third  central  moment;  can  be  used  in 
testing  for  normality 

BETA2:   Biased  estimate  of  fourth  central  moment 

MAXIMUM:   Largest  sample  value 

MINIMUM:   Smallest  sample  value 

QUANTILES:   The  al pha-quanti 1 e,  x(alpha),  is  the  solution  to  the 
equation: 

Probability  (x  .le.  x(alpha)  )  =  alpha 
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ship  test  days;  and  4)  the  GE0CEI7ER  STl  31370  position  and 
the  uncorrected  stationary  GPS  position  for  the  three  days 
of  ship  tests.  In  addition,  the  azimuths  between  the 
theodolite-computed  static  corrections,  for  each  subset, 
were  graphed  as  a  series  of  histograms.  To  illustrate  the 
difference  between  the  GE0CEI7ER  STA  31370  position  and  the 
stationary  GPS  position,  the  azimuths  between  the  two 
positions  are  graphed  on  histograms  for  each  of  the  three 
test  days. 
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71.   ERROR  ANALYSIS 
The  error  sources  contributed  by  each  positioning  system 
are  described.   Estimates  of  the  effect  of  each  on  the  final 
results  are  determined  where  possible. 

A.   ARTEMIS  POSITIONING  SYSTEM 

The  sources  of  error  in  the  ARTEMIS  ship  position  can  be 
divided  into  two  categories — those  sources  of  error  inherent 
in  all  microwave  positioning  systems  and  those  unique  to  the 
ARTEMIS  system.  All  are  investigated  and  the  contribution 
of  each  to  the  total  error  is  determined. 
1 ,   Microwave  System  Errors 

Sources  of  system  error  include  range  holes, 
inaccuracies  in  calibration,  non-ideal  geometry  of  the  lines 
of  position  describing  the  fix,  and  interference  with 
reflections  of  the  signal  from  nearby  objects. 

Range  holes  are  a  cause  of  unreliability  in 
microwave  positioning  systems.  The  problem  exists  when  the 
sea  surface  is  calm  enough  to  reflect  part  of  the 
transmitted  signal.  The  reflected  signal  will  arrive  at  the 
receiving  antenna  at  a  later  time  than  the  signal  which 
travels  directly,  causing  a  phase  difference  between  the  two 
signals.    A  zone   or  region  of  destructive   interference  of 
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the  directly  transmitted  signal  and  the  reflected  signal 
will  occur  within  which  the  vessel  will  receive  an 
attenuated  signal.  When  the  phase  difference  of  the  two 
paths  is  exactly  180  °,  no  signal  will  be  received  provided 
the  reflected  signal  is  as  strong  as  the  direct-path  signal. 
Prediction  of  the  range  holes,  therefore,  is  useful  in 
survey  planning  so  as  to  avoid  this  problem. 

The  range  hole  is  a  function  of  shore  unit  and  ship 
unit  antenna  heights,  the  wavelength  of  the  system,  and  the 
distance  from  the  shorebased  unit.  If  the  operating  area  is 
known,  antenna  heights  can  be  determined  and  equipment  can 
be  installed  so  that  no  range  holes  exist  in  the  survey 
area.  If  fixed  antenna  heights  are  known,  the  region  of 
uncertainty  can  be  determined  from  the  equation:  [Bef.  48] 
The  following  equation  is  used  for  this  computation: 

Dp  =  2(h1)  (h2)/nX, 

where  X  is  the  wavelength  (constant  for  ARTEMIS  at  .032  m) , 
hi  and  h2  are  antenna  heights,  and  D^  is  the  distance  from 
the  shore  based  unit  to  the  nth  range  hole.  For  this  study, 
the  antenna  height  of  the  fixed  unit  above  mean  sea  level  at 
the  GE0CSI7ER  STA  31370  was  3.8  meters.  The  mobile  antenna 
was  installed   on  top  of  the   equipment  shelter  on   the  boat 
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deck  of  the  R/V  &CANIA  at  a  height  of  7.9  meters  above  the 
the  waterline.  Two  range  holes  within  the  survey  area  were 
determined  to  be  located  at  ranges  of  1376  meters  for  n  =  1 
and  938  meters  at  n  =  2  (fig.  6.1).  These  two  ranges  were 
computed  for  the  case  of  mean  lower  low  tide.  The  change  in 
width  of  the  range  hole  due  to  a  tidal  fluctuation  of  1.5 
meters  during  the  test  is  shown  in  figure  6,1.  (During  the 
night  of  15  May  the  sea  surface  was  sufficiently  disturbed 
by  a  two-foot  swell  and  a  surface  wind  of  2-4  knots  such 
that  surface  reflection  was  minimal  and  no  loss  of  signal 
was  observed  at  the  computed  ranges) . 

Calibration  of  the  ARTEMIS  system  was  accomplished 
in  two  stages.  As  part  of  the  study,  the  ARTEMIS  was 
obtained  for  a  one  week  test  of  the  system  from  13  to  17 
April.  During  this  time  it  was  calibrated  for  azimuth  and 
range  over  the  baseline  between  the  GE0CEI7EH  STA  31370  and 
MOSSEL.  The  azimuth  readings  required  no  adjustment  but  the 
range  was  reduced  by  two  meters.  The  equipment  was  returned 
to  Marinav  Corporation  to  fulfill  prior  commitments  and  then 
was  returned  to  Monterey  for  the  ship  test  on  15  May. 
Before  the  evening  operation,  the  ARTEMIS  was  used  to 
measure  the  same  baseline  as   previously.    A  plus  two  meter 
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bias  was  observed  at  this  time,  but  no  immediate 
compensation  was  made.  This  error  was  assumed  constant 
throughout  the  test  and  was  included  in  the  calculation  of 
the  ship-determined  ARTEMIS  position  during  postprocessing. 

One  of  the  advantages  of  the  ARTEMIS  system,  unlike 
other  microwave  positioning  systems,  is  that  the  geometry  is 
always  ideal  for  the  lines  of  position  that  determine  the 
fix.  The  range-azimuth  technique  of  determining  the  vessel 
position  defines  an  uncertainty  ellipse  having  axes  that  are 
always  mutually  perpendicular.  This  allows  prediction  of 
errors  as  a  function  of  range  only,  independent  of  the 
angle. 

The  manufacturer's  specified  error  in  range 
measurement  of  the  ARTEMIS  is  ±  1.5  meters  for  all  line  of 
sight  ranges.  One  axis  of  the  uncartainty  ellipse  is 
therefore  constant.  The  other  axis,  that  describes  the 
error  introduced  by  the  system  azimuth  uncertainty,  is  60  cm 
per  1000  m  range.  (This  is  determined  from  an  angular  error 
of  0.033°  or  2'  of  arc).  Both  system  error  values  represent 
an  accuracy  of  twice  the  standard  deviation  value,  or  a  95'5 
probability  of  the  true  value  occurring  within  the  margin 
allowed  by  these  parameters. 
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The  positional  uncertainty  can  be  expressed  in  terms 
of  the  radius  of  a  circle  of  equivalent  probability  within 
which  the  true  fix  is  expected  to  occur.  The  method  of 
computation  of  this  radius  is  described  in  Appendix  A.  The 
technique  allows  the  conversion  of  uncertainty  ellipses  to 
radii  of  circles  of  equivalent  probability.  Values  from 
this  computation  for  the  radii  of  circles  of  90%  probability 
at  different  ranges  used  in  the  survey  area  for  15  May  are 
displayed  in  figure  6.2  along  with  the  circlss  of  90% 
probability,  the  magnitude  of  the  radii  at  the  ranges  used 
for  the  study  do  not  exceed  5.4  meters.  One  goal  of  this 
study  was  to  measure  the  difference  in  positional  accuracy 
between  ARTEMIS  and  GPS.  A  maximum  error  of  1.0  meters,  an 
order  of  magnitude  smaller  that  the  tested  system  accuracy, 
is  desirable.  At  the  present  time,  however,  there  are  few 
systems  that  can  provide  the  order  of  accuracy  system  for 
the  position  of  a  moving  vessel,  so  the  5.4  meter  error 
observed  at  4000  meters  had  to  suffice. 

Error  caused  by  reflection  of  the  signal  from  nearby 
objects  is  difficult  to  measure.  To  minimize  this  problem, 
the  offset  position  for  the  ARTEMIS  unit  was  planned  so  as 
to   be  clear   of  possible   interference   from   the   wooden 
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structure  and  the  ej uipment  vehicle.  No  signal  interference 
was  observed  despite  the  ship's  motion.  This  source  of  data 
unreliability  was  therefore  not  considered  to  be 
significant. 

2.   ARTEMIS  System  Errors 

Errors  unique  to  the  ARTEMIS  system  include  a 
pointing  error  source  and  the  lack  of  resolution  of  the 
azimuth  measurement  in  the  dynamic  mode  of  operation. 

The  operation  of  the  ARTEMIS  system  requires  a 
manual  alignment  of  the  fixed  antenna  so  that  it  is 
perpendicular  to  a  known  object.  This  is  necessary  to 
establish  an  original  reference  direction  upon  which  all 
angular  measurements  are  based.  It  is  possible  that  a 
pointing  error  could  be  produced  by  the  operator  if  the 
crosshairs  of  the  telescope  used  to  align  the  antenna  are 
not  exactly  sighted.  This  would  introduce  a  systematic 
error  into  all  subsequent  angle  measurements.  To  minimize 
this  source  of  error,  the  alignment  of  the  antenna  was 
checked  before  and  after  the  period  of  operation.  The 
original  azimuth  displayed  on  the  antenna  unit  display 
agreed  exactly  with  the  final  azimuth  value.  Any 
contribution  to  ship  position  error  was  considered  minimal. 
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The  AHTEHIS  was  operated  using  the  dynamic  mode 
option  (see  Squipment  Descript ion--ARTEMIS) .  The  resolution 
of  the  system  in  this  mode  is  1  meter  for  displayed  range 
and  0.1°  for  azimuth  values.  This  resolution  limit 
introduces  an  uncertainty  of  _0.5  meters  m  range  and  0.05 
in  angle  measurements.  The  computed  uncertainty  ellipse  as 
expressed  by  radii  of  circles  of  equivalent  probability,  is 
shown  in  figure  6.2.  The  system  error  and  all  other  errors 
mentioned  above  are  assumed  to  occur  within  the  scope  of 
this  resolution  uncertainty. 

Another  source  of  system  error  has  been  proposed  in 
a  study  conducted  in  June,  1980,  by  NOAA.  [Ref.  49]  The 
test  compared  the  ARTEMIS  system  with  two  ranges  measured  by 
two  Del  Norte  Trisponders  (a  short  range  microwave 
range/range  positioning  system)  and  two  lines  of  direction 
measured  by  two  shore-based  theodolites.  The  results 
described  an  underway  bias  of  -0.05°  of  unknown  origin  in 
the  angle  measurements.  This  discrepancy  has  aot  been 
resolved.  The  preliminary  test  conducted  in  Monterey  in 
April,  1981  compared  the  ARTEMIS  system  with  three  lines  of 
direction  determined  by  three  theodolites.  The  purpose  of 
the  test  was  to  resolve  the  underway  bias  problem.  Results 
are  enclosed  in  a  separate  report.   [Ref.  50] 

in 


B.   CONTROL  NETWORK 

The  entire  survey  was  oriented  relative  to  the  TRANSIT 
position  on  the  Coast  Guard  Pier  (GE0C2I72R  STA  31370). 
Therefore,  errors  inherent  in  the  position  of  this  station 
mark  affected  all  results.  Assuming  that  the  geographic 
coordinates  published  by  DSA/HTC  represented  the  best 
estimate  available  of  the  GEOCEIVER  STA  31370,  the  remaining 
significant  sources  of  error  affecting  all  results  were 
survey  errors.  The  control  network  for  the  study  was 
determined  from  the  traverse  and  it  was  assumed  that  all 
positions  used  in  the  study  were  oriented  relative  to  each 
other  on  the  HGS-72  ellipsoid.  A  change  in  the  ellipsoidal 
surface  within  the  area  encompassed  by  the  traverse  was 
assumed  to  be  insignificant. 

Survey  errors  stem  from  the  methods  used  for  obtaining 
the  initial  azimuth  for  orientation  of  the  survey,  the  type 
of  traverse  chosen,  and  the  geometry  of  the  traversed 
triangle. 

To  obtain  a  starting  azimuth  for  the  traverse,  two  sets 
of  observations  on  Polaris  were  recorded  at  the  GEOCEIVER 
STA  31370.  A  striding  level  attachment  to  the  theodolite 
was  not   used.    The   instrument  level   bubble  was   observed 


112 


repeatedly  during  the  observations  to  ensure  constant 
instrument  alignment  with  the  vertical.  Failure  to  keep  the 
instrument  level  could  have  introduced  a  +  30"  uncertainty 
into  the  starting  azimuth  of  the  survey.  [Ref.  51]  To 
consider  the  significance  of  such  an  error,  all  stations 
were  recomputed  using  a  deviation  of  +  15"  from  the  computed 
azimuth.  The  resulting  shift  in  the  positions  of  MUSSEL  and 
aSE  MON  are  listed: 

USE  MON      0.133  m  distance  at  230° 36*47 .  16"  azimuth 

MUSSEL      0.132  m  distance  at  230° 33*  4.08"  azimuth 

The  change  is  considered  insignificant  relative  to  the  other 
errors  described  above. 

The  survey  was  a  closed  loop  traverse  using  3E0CEIVSR 
STA  31370  as  the  origin  and  closing  position.  The  weakness 
in  this  method  involves  the  measured  distances.  A  method  of 
checking  the  observed  angles  in  the  triangle  exists  via  the 
angular  closure  if  a  significant  deviation  in  one  of  the 
angles  can  be  detected.  The  distance  measurement,  however, 
can  produce  an  undetected  bias  in  the  survey  if  the  distance 
measuring  instrument  contains  a  systematic  error.  This 
error   would  shift   all   the   surveyed  positions.     Correct 
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measurement  of  distance  by  the  MRA  V  Tellurometer, 
therefore,  was  important.  The  equipment  was  calibrated  over 
a  known  baseline  and  the  errors  were  within  the  centimeter 
range.   This  error  was  therefore  negligible. 

The  geometry  of  the  triangle  was  not  optimum.  Ideally, 
the  triangle  would  have  been  equilateral  to  minimize  the 
source  of  observing  error.  The  choice  of  traversed  points 
was  based  on  the  location  of  preexisting  station  marks  and 
the  desired  geometry  of  the  theodolite  fixes  to  the  ship. 
The  closures  computed  for  the  traverse  met  third  order 
specifications  for  accuracy,  and  the  survey  was  not  adjusted 
to  remove  the  closure  error.   [Ref.  52,  53J 

The  survey  errors  were  considered  to  be  within  the 
uncertainty  of  the  TRANSIT  station  location  and  were 
therefore  not  considered  significant  for  the  study. 

C.   THEODOLITE  NETWORK 

The  errors  generated  by  the  theodolite  network  are 
divided  into  equipment  errors,  operator  related  errors  and 
system  errors,  of  which  the  latter  is  the  most  significant. 

The  equipment  used  for  the  theodolite  network  were  three 
1-second  theodolites  which  have  a  measurement  resolution  of 
10.5".     All   theodolites    were   in   good   condition. 
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Instrumental  errors  of  a  theodolite  are  listed  in  Table  VI. 
[Ref.  54] 

The  lajor  source  of  equipment  error  was  introduced  by 
not  changing  face,  the  process  of  observing  angles  a  second 
time  Mith  the  telescope  reversed.  Due  to  the  nature  of  the 
test  coupled  with  1-minute  fix  intervals,  it  was  not 
possible  to  observe  all  the  fixes  in  both  the  direct  and 
reversed  modes.  As  an  example  of  the  magnitude  of  this 
error,  the  horizontal  collimation  errors  were  determined  for 
each  theodolite  prior  to  the  study.  The  results  are  listed 
below: 

Theodolite  Collimation  5rroip 

»  iauo5  -  1.0" 

#  14482  *    2.7" 

»  14452  -  0.7" 

All  the  measured  angles  were  adjusted  for  this  error  as 
described  in  the  processing  section.  The  largest  influence 
on  the  remaining  instrumental  errors  was  the  non-vert icality 
of  the  main  axis.  If  the  non-verticality  were  large  (for 
example,  one  division  off-level  in  the  lower  spirit  level 
bubble) ,  the  angle  measured  could  be  in  error  by  7". 
[Ref.    55]        To    minimize     this      error,        each    theodolite      was 
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TABLE  VI.  Theodolite  Instrument  Errors 


ERROR 


ELIMINATED  BY 


Eccentricity,  misplaced  indices 

Graduation  of  circles 
Micrometer  run  and  vernier  error 

Non-verticality  of  main  axis 


Bad  setting  of  the  upper  spirit 
level 

Vertical ity  of  cross  wires 


Horizontal   collimation 

Vertical  collimation  and  index 
setting  error 


Reading  both  sides  of  circle 
and  changing  face 

Changing  zeros 

Selection  of  suitable  zero 
setting 

Error  is  not  cancelled  by  any 
observing  procedure 

Correct  setting  or  applica- 
tion of  correction 

Changing  face  and  using  same 
part  of  wire 

Changing  face 

Changing  face 
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relevelled  if  necessary  between  sets  of  observations.  The 
magnitude  of  this  and  the  collimation  errors  were  considered 
to  be  insignificant  relative  to  the  system  errors  discussed 
later  in  this  section. 

The  operator-related  errors  are  a  function  of  the 
individual  performance  of  the  observers.  The  experience  and 
personal  comfort  of  the  observer,  for  example,  are  two  major 
factors  which  can  influence  the  observed  angles.  An 
observer's  performance  and  resultant  impact  on  the  guality 
of  the  results  depend  on  the  judgment  and  reaction  time  of 
the  observer  and  are  influenced  by  variables  such  as  fatigue 
and  age.  The  subjective  nature  of  the  problem  makes  an 
estimate  of  the  error  difficult,  therefore,  for  this  study 
the  error  was  considered  random. 

The  major  source  of  error  in  the  theodolite  positioning 
system  dealt  with  the  timing  of  the  fix.  A  delay  existed 
between  the  time  the  • mark*  was  called  over  the  radio  and 
the  actual  time  the  angle  was  observed.  The  timing  problem 
was  most  critical  when  the  ship  was  close  to  one  of  the 
three  observer  locations  because  the  rate  of  change  of  the 
bearing  of  the  ship  was  greatest  at  this  time.  Hence  a 
small  delay  generated  a  large  amount  of  error  in  the  angle 


117 


recorded.   To  indicate  the  effect  this  error  source  may  have 
had  on  the  position   of  the   ship,   two   fixes  were   chosen 
arbitrarily — one  at   close  range  to  each  of  MUSSEL   and 
GEOCEIVEH  STA  31370.     A  one  ninute  value  was   added  to  the 
angle  which   was  most  affected  by  the  timing  error   and  the 
fixes  were  recomputed.     The  same  process  was  repeated  for 
the  three   fixes  with  the  one  minute  value   subtracted  from 
the  critical  angle  in  each  fix.   The  results  follow: 
Number    Original  Fix    *•  1  Minute       -  1  Minute 
210       36 °37« 20.62282"   36° 37» 20.6 1468"   36°37» 20. 63081 " 
121°53«iH.  9718U"  121°  53»  4  1 .96933"  1  21  °53' U  1 .  97U26" 
259       36°36»45. 88048"   36°  36« 45.87574"   36 °36 • 45. 89082" 
121°52«54. 05062"  121  °  52» 54. 04635"  121°52« 54. 05459" 
No  significant  error  was  generated  from  this  source. 

Another  source  of  system  error  was  associated  with  the 
geometry  of  the  fix.  Ideally,  the  three  theodolite  angles 
would  describe  two  equilateral  triangles  to  the  ship; 
however,  within  the  operating  area,  the  geometry  varied.  As 
an  indication  of  the  size  of  this  error,  representative 
fixes  were  chosen  and  the  circular  map  accuracy  standards 
(CMAS)  of  each  of  the  combinations  was  determined.  Figure 
6.3  illustrates  the  approximate  error  caused   by  variations 
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in  the  geometry.  Within  the  area  outlined  in  the  figure  the 
geometry  creates  a  CMAS  of  not  more  than  10. U  meters.  This 
is  therefore  the  source  of  the  largest  probable  error  in  the 
theodolite  method  of  positioning.  Nevertheless,  this  method 
in  combination  with  the  ARTEMIS  system  provided  the  best 
estimate  of  the  ship  location. 

D.   TRANSIT  POSITIONING  STSTEH 

The  method  used  to  determine  the  TRANSIT  position 
(GE0CEI7EE  STA  31370)  is  the  most  accurate  satellite 
positioning  technique  presently  available.  Error  sources 
are  reduced  both  during  data  collection  and  during 
postprocessing  of  the  recorded  values. 

DMA  places  requirements  on  the  collection  of  data  that 
will  minimize  error.  For  this  study,  only  satellite  passes 
whose  maximum  altitude  ranged  between  10°  and  30  ^  were 
recorded.  Meteorological  data  (temperature,  pressure,  and 
humidity)  were  collected  for  each  pass,  and  a  determination 
of  the  receiver  clock  drift  was  computed  after  each  pass. 
The  amount  of  clock  drift  as  a  function  of  time  gave  a 
real-time  indication  of  receiver  performance.  Equipment 
checks  were  conducted  routinely,  and  the  results  were 
included  with  the  data  package  submitted  to  DMA. 
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An  error  budget  of  the  TRANSIT  system  for  point 
positioning  is  listed  in  Table  VII.  [Ref.  56]  The  precise 
ephemeris  method  of  postprocessing  accounts  for  all  of  these 
errors:  1)  for  each  pass  the  tropospheric  correction  is 
computed  from  the  meteorological  data;  2)  the  ionospheric 
correction  is  determined  by  observing  the  shift  in  frequency 
between  the  two  transmitted  frequencies;  3)  the  receiver 
clock  error  is  known;  4)  the  satellite  orbits  are  observed 
by  the  tracking  stations  so  that  precise  satellite  positions 
are  known;  5)  errors  due  to  uncertainties  in  the 
geopotential  model  and  surface  drag  forces  are  accounted  for 
during  computation;  and  6)  the  receiver  altitude  relative  to 
the  ellipsoidal  surface  is  computed  approximately  from  the 
given  MSL  elevation. 

The  accuracy  of  the  computed  TEANSIT  position  is  stated 
as  1.5  meters  per  axis  within  a  90%  confidence  interval. 
This  is  assumed  to  provide  the  best  estimate  possible  of  the 
GEOCSIVER  STA  31370  relative  to  the  MGS-72  coordinate 
system. 

E.   GPS  ERROR  SOURCES 

Global  Positioning  System  error  sources  are  separated 
into  two  main   categories,   the  space  vehicle   systam  errors 
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TABLE  VI I. TRANSIT  System  Errors 


SOURCE  ERROR  (m) 

Uncorrected  propagation  effects 
(ionospheric  and  tropospheric  effects)       1-5 

Instrumentation  and  measurement  noise 
(local  and  satellite  oscillator  phase 
jitter,  navigator's  clock  error)  3-6 

Uncertainties  in  the  geopotential  model 
used  in  generating  the  orbit  10-20 

Incorrectly  modelled  surface  forces 
(drag  and  radiation  pressure  acting  on 
the  satellites  during  the  extrapolation 
interval)  10-25 

Uncertainties   in  navigator's  altitude 
(results   in  a  bias   in  longitude)  10 

Ephemeris   rounding  error  5 
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and  the  GPS  user  receiver  errors.  There  are  seven  possible 
sources  of  error  in  the  space  vehicle  systein--errors 
associated  with  satellite  equipment,  satellite  clock, 
ephemeris,  atmospheric  delays,  signal  multipath  effect, 
satellite  update,  and  satellite  propulsion  firing.  Table 
III  lists  these  error  sources  and  the  approximate  user 
equivalent  range  error  for  each.  [Bef.  57] 

The  satellite  equipment  error,  also  termed  the  group 
delay,  is  caused  by  differences  and  uncertainties  in 
processing  and  transmission  of  the  message  through  the 
individual  satellite's  circuitry.  Each  satellite  has  the 
group  delay  calibrated  prior  to  launch  and  this  known  group 
delay  is  included  in  the  satellite  time  corrector  within  the 
navigation  message.  However,  unknown  delays  in  signal 
processing  and  circuitry  not  accounted  for  by  the  corrector 
contributes  an  estimated  one  meter  to  the  error  budget. 

Satellite  clock  error  is  caused  by  variations  in  each 
satellite  clock.  Satellite  clock  time  may  depart  by  up  to 
976  microseconds  from  the  correct  GPS  time.  Corrections  for 
the  predicted  clock  variations  are  included  in  the 
navigation  message.  Residual  clock  variations  contribute  to 
the  error  budget  and  are   included  in  ephemeris  errors  since 
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the  clock  errors  are  small  and  have  the  same  affect  as 
ephemeris  errors. 

Ephemeris  errors  are  the  result  of  satellite  clock 
variations,  drift  by  the  monitor  station  clocks,  monitor 
station  signal  processing  delays,  lack  of  precise  monitor 
station  geodetic  positions,  solar  wind  and  pressure 
parameter  variations,  and  earth  gravity  model  error.  These 
sources  all  contribute  to  the  ephemeris  error  budget.  The 
ephemeris  error  budget  is  monitored  over  long  time  periods 
and  translated  by  the  Master  Control  Station  into  an 
Oser-Eguivalent  Range  Error  (ORE).  When  the  ORE  exceeds  the 
prescribed  value  of  4  meters  the  satellite  ephemeris  data  is 
uploaded  to  reduce  the  ORE.  [Bef.  58]  Most  of  the  ORE  is 
included  into  the  user  clock  corrector  obtained  in  the  GPS 
position  computation.  However,  not  all  of  the  ephemeris 
error  is  removed  during  the  position  computation,  and  the 
residual  error  has  been  estimated  to  be  approximately  1.5 
meters.   See  Table  III. 

Atmospheric  delays  in  signal  propagation  are  the  result 
of  two  effects,  ionospheric  delay  and  tropospheric  delay. 
Of  these,  the  ionospheric  delay,  which  depends  on  the 
freguency  of  the  transmitted  signal,   is  greater.    With  the 
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two  frequencies,  LI  and  L2,  transmitted  by  each  satellite  to 
the  user  receiver,  the  ionospheric  propagation  delay  can  be 
modelled.  Osing  this  model,  the  majority  of  the  propagation 
delay  is  accounted  for,  but  a  residual  delay  results  in  an 
error.  The  other  part  of  the  atmospheric  delay  error 
depends  on  the  tropospheric  propagation  delay.  Propagation 
delay  in  the  troposphere  is  a  function  of  signal  path  length 
through  the  layer,  the  receiver  air  temperature,  and  air 
pressure  and  humidity.  Since  the  tropospheric  delay  is  not 
a  function  of  frequency  and  is  highly  localized,  the 
tropospheric  delay  corrector  has  not  been  incorporated  into 
the  atmospheric  correction  computation  and  is  therefore  an 
additional  source  of  error.  These  two  error  sources  combine 
to  yield  an  atmospheric  delay  error  estimated  to  have  a 
magnitude  of  between  2.U  to  5.2  meters. 

Ionospheric  scintillation  of  the  GPS  signal,  in  addition 
to  the  ionospheric  propagation  delay,  cause  errors  in 
positioning.  Scintillation  has  the  effect  of  random 
interference  with  the  propagated  signal  as  the  signal  is 
transmitted  through  the  upper  ionosphere.  [Hef.  59] 
Interference  caused  by  scintillation  combines  with  the  GPS 
user  receiver   processing  noise  to  cause  a  random   ranging 
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error  and,  if  the  interference  is  severe,  interruption  of 
the  pseudo-random  number  binary  code  causes  loss  of  signal. 
Errors  and  signal  losses  caused  by  scintillation  cannot  be 
predicted  or  estimated. 

Signal  nultipath  errors  result  from  the  reception  of  a 
satellite's  signal  transmitted  via  different  propagation 
paths.  Received  signals  thus  may  have  been  distorted  or 
reduced  in  strength  to  a  sufficient  degree  so  as  to  affect 
the  range  determination  from  the  satellite.  Error  due  to 
multipathing  is  highly  localized  and  cannot  be  modelled. 
Estimates  of  the  error  magnitude  are  between  1.2  and  2.7 
meters. 

The  last  two  space  vehicle  system  error  sources  are 
attributable  to  the  Master  Control  Station  (MCS)  satellite 
accuracy  and  maintenance  functions.  Usually  the  MCS  uploads 
each  satellite  at  a  given  time  each  day  and  the  ephemeris 
data  is  changed(Table  7III) .  Should  the  user  not  be  aware  of 
occurrence  of  upload  and  if  the  receiver  therefore  does  not 
obtain  the  new  ephemeris  data,  ranging  to  the  satellite  in 
question  will  contain  a  possibly  very  large  error.  In 
addition,  if  the  HCS,  through  the  Remote  Stations, 
determines   that  the   four  skewed   reaction   wheels   cannot 
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maintain  satellite  antenna  alignment  with  earth,  the 
hydrazine  propulsion  system  is  activated  and  the  satellite 
is  rotated.  Once  the  propulsion  system  has  used,  the 
satellite  orbit  cannot  stabilize  immediately,  and  the 
ephemeris  data  will  be  inaccurate.  A  period  of  up  to 
several  hours  may  be  required  for  orbit  stabilization  and 
for  the  new  ephemeris  data  to  be  computed  and  uploaded. 
During  that  period,  the  ranging  computations  to  the 
satellite  will  have  large  errors.  The  magnitude  of  these 
two  error  sources  cannot  be  estimated  but  is  assumed  to  be 
several  orders  of  magnitude  larger  than  the  previously 
stated  errors. 

The  second  error  category  is  that  related  to  the  GPS 
user  receivers.  These  sources  of  error  associated  with  the 
user  GPS  receiver,  the  Manpack,  are  receiver  signal  noise, 
receiver  vessel  velocity,  signal  resolution,  method  of 
signal  processing,  ephemeris  update  capability,  and,  for  the 
differential  GPS  computations,  system  accuracy  comparison. 

Received  signal  noise  errors  result  from  interference 
introduced  into  the  receiver  GPS  signals  through  processing 
and  internal  receiver  circuitry.  Additionally,  quality  of 
signal  resolution  is  dependent  upon  receiver  type  and  can  be 
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a   source  of   error.    The  combined  errors   from  these   two 
sources   are  estimated  to   be   1.5   meters.    [Ref.  60] 
Estimates  of  these  error  sources  for  the  avUE  have  not  been 
made  but  are  estimated  to  be  of  the  same  magnitude. 

Velocity  of  the  receiver  vessel  results  in  GPS  errors  of 
varying  magnitude.  The  higher  the  velocity,  the  larger  the 
error.  This  error  is  due  to  the  time  lag  required  to  obtain 
and  process  each  NAVSTAR  satellite  range  with  the  MVOE 
single  channel  receiver,  to  compute  the  average  vessel 
velocity  and  to  use  the  average  velocity  for  positional 
computations.  The  Manpack  requires  a  six  second  time 
interval  to  obtain  and  process  the  pseudo-range  from  each 
satellite  used.  Thus  the  fix  position  is  an  extrapolation 
of  the  four  satellite  ranges  to  obtain  a  convergent  fix 
position.  At  one  knot,  the  vessel  travels  0.51  meters  in 
one  second.  During  the  full  twenty-four  seconds  required  to 
obtain  four  satellite  ranges,  the  vessel  would  have 
travelled  12. 2U  meters.  With  the  vessel  travelling  at  four 
knots,  the  distance  travelled  during  the  total  time  interval 
required  for  a  GPS  fix  position  would  be  48.96  meters. 
Processing  the  pseudo-ranges  from  four  satellites  requires 
the  solution   including  the   average  vessel   velocity  during 
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the  2U-second  ranging  period.  Use  of  this  average  neglects 
the  relative  velocity  of  the  recaiver  antenna  which  is 
mounted  on  the  mast.  During  this  test,  the  GPS  antenna  was 
mounted  9.06  meters  above  the  waterline  of  the  R/V  ACANIA. 
Assuming  five  degrees  of  pitch  and  roll  and  a  vessel  roll 
period  of  15  seconds,  relative  average  antenna  velocity  was 
,0.16  meters  per  second.  The  antenna  velocity  error 
translates  into  a  random  error  of  0.16  meters.  The  other 
velocity  errors  have  not  been  estimated  but  are  assumed  to 
be  less  than  one  meter  for  the  low  vessel  velocities  used  in 
the  test. 

One  error  source  mentioned  previously  requires  further 
clarification:  error  due  to  the  relative  position  at  the 
time  of  the  fix.  The  ManpacJc  requires  up  to  twenty-four 
seconds  to  obtain  and  process  the  satellite  ranges.  The 
position  of  the  fix  and  the  associated  time  of  fix  are  then 
recorded.  However,  the  time  of  the  fix  may  not  have  been 
the  same  time  the  computational  results  were  printed.  Thus 
the  output  position  may  not  have  bsen  correlated  with  the 
correct  theodolite  computed  ship  position  since  the 
theodolite  lines  of  position  were  observed  when  the  GPS 
position  was   printed.    Errors  from  this   bias  source  may 


130 


range  from  12.24  meters  at  one  knot  of  vessel  velocity  to 
43.96  meters  at  four  knots.  This  error  varied  with  the 
processing  time  required  to  obtain  a  convergent  position. 
An  estimate  of  this  error  magnitude  was  not  made  because  of 
its  wide  variability. 

The  type  of  receiver  and  the  method  of  signal  processing 
also  contributed  to  the  error.  Multichannel  receivers 
designed  for  Phase  I  testing  have  a  lower  magnitude  error 
for  satellite  signal  reception  and  processing  than  do  the 
single  channel  receivers,  such  as  the  Manpack.  In  the 
multichannel  receivers,  four  or  mors  satellite  pseudo-range 
data  are  taken  simultaneously.  The  single  channel  receiver 
sequentially  receives  each  satellite  and  computes  the 
pseudo-range  data.  Simultaneous  data  reception  reduces  the 
user  velocity  error  and  therefore  has  less  position  aliasing 
than  that  found  in  the  single  channel  receivers. 

The  method  of  signal  processing  may  induce  positional 
errors  through  the  use  of  non-optimal  weighting  functions  in 
the  position  computation.  Should  the  ephemeris  or  ranging 
data  change  significantly,  as  happens  during  uploading, 
signal  processing  using  Kalman  filters  as  used  in  the 
aanpack  do  not  change  the  computed  position   rapidly  enough 


131 


to  reflect  the  large  changes  in  ranges  and  result  in  a  slow 
change  drift  of  the  GPS  position.  Estimated  time 
requirements  for  user  receivers  to  stabilize  after  an  upload 
are  up  to  two  hours.   [Ref.  61] 

The  capability  of  the  user  GPS  receiver  to  receive 
ephemeris  data  immediately  after  upload  reduces  the 
ephemeris  errors  caused  by  using  an  outdated  ephemeris.  The 
ephemeris  data  is  repeated  every  thirty  seconds  and  its 
length  requires  sections  of  two  message  subframes. 
( Bef  •  62,  630.  The  use  of  current  ephemeris  data  in 
the  computations  reduces  the  error  for  the  ephemeris  to  the 
value  shown  in  table  III.  The  Manpack  has  the  capability  to 
accept  new  ephemeris  data  immediately,  but  during  the  test 
large  errors  were  introduced  into  the  computation  on 
satellite  ephemeris  upload  because  this  selected  optional 
function  either  did  not  operate  properly  or  was  neglected  by 
the  operators. 

?or  the  differential  GPS  user,  an  additional  error 
source  is  introduced  into  the  error  budget:  the  assumption 
that  both  GPS  receivers  will  observe,  obtain  the 
pseudo- ranges ,  and  compute  the  position  solution  in  exactly 
the  same   manner.    If   this  were   true  then   both  receivers 
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should  be  able  to  compute  the  identical  position  if  placed 
at  the  same  position  at  the  sane  time.  With  the  Manpacks, 
the  assumption  that  the  units  were  identical  and  the 
position  solutions  were  identical  was  not  made  because  prior 
testing  established  that  the  receivers  differed  in  position 
computations  by  as  much  as  3.08  meters.   [Ref.  64] 

For  the  static  testing  on  11  May,  one  error  source  not 
previously  discussed  was  the  vertical  separation  of  the  two 
Manpaclc  antennas,  colocated  over  the  GEOCEIVER  ST  A  31370 
control  station,  and  the  possible  interference  of  the 
uppermost  antenna.  System  1,  with  the  GPS  signals  received 
by  the  lower  System  2  antenna.  Separation  of  the  two 
antennas  was  2.1  meters.  Calculations  using  the  speed  of 
propagation  in  a  vacuum  and  the  two  GPS  frequencies  produce 
wave  lengths  of  0.19  meters  and  0.24  meters.  Since 
separation  of  the  antennas  was  much  greater  than  either 
wavelength  theoretically  there  should  have  been  no 
positional  errors  due  to  antenna  interaction  or 
interference.  Differing  elevations  of  the  antennas  above 
the  ellipsoidal  surface  also  may  have  induced  errors  in  the 
position  computation.  Estimated  errors  for  the  elevation 
differences  are  less  than  0.1  meters.    Total  errors  for  the 
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antenna   vertical   elevation    separation   and   antenna 
interference  or  interaction  have  not  been  estiaated. 

Thus,  the  error  budget  for  the  Global  Positioning  System 
suggests  a  combined  error  of  up  to  57.35  meters.  tJsing  the 
differential  technique,  the  bias  errors  can  be  reduced  in 
magnitude  or  eliminated.  These  bias  errors  include 
ionospheric  propagation  delays,  tropospheric  propagation 
delays,  variations  from  predicted  ephemerides,  and  satellite 
clock  perturbations.  Teasley,  Hoover,  and  Johnson  (1980), 
evaluated  the  error  budget  for  a  more  advanced  GPS  receiver 
and  found  the  predicted  filtered  solution  error  was  reduced 
from  approximately  11  meters  to  approximately  3  meters  with 
differential  GPS  processing.  [Ref.  65]  Error  budgets  for 
the  aanpaclc  and  the  differential  correction  vector 
processing  method  used  in  this  test  did  not  attain  the  four 
meter  accuracy  but  were  not  expected  to  because  the 
receivers  used  were  less  accurate,  and  of  different  types. 
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VII.  RESULTS 
Results  of  the  GPS  tests  have  been  separated  into  four 
sections:  stationary  GPS  evaluation,  evaluation  of  non- 
differential  GPS,  differential  GPS  evaluation  with  static 
correction,  and  evaluation  of  diffsrential  GPS  without 
static  correction.  Tables  IV  and  VIII  list  the  optimum 
satellite  observation  periods  and  the  NAVSTAR  upload  times 
for  each  day  of  the  testing  period, 

A.   STATIONARY  GPS  EVALUATION 

Stationary  GPS  receiver  evaluation  was  separated  into 
the  stability  tests,  the  static  tests  and  the  tests  of  the 
stationary  GPS  receiver  used  during  differential  testing. 
Figures  7.1,  7.2,  and  7.3  illustrate  the  results  for  the 
stability  test  conducted  on  9  and  10  May.  The  receiver  was 
less  stable  than  expected  and  depended  to  a  great  extent  on 
the  HAVSTAR  uploading  times  and  the  UHE  of  NAVSTAR  1  (fig. 
7.4  and  7.5) .  Correlation  between  times  for  the  NAVSTAR 
upload  and  changes  in  the  stability  can  be  observed  at  0620 
on  9  May  (fig.  7.1),  0655  on  9  May  (fig.  7.2),  0815  on  10 
May  (fig.  7.3),  and  0845  on  10  May  (fig.  7.3).  The 
relationship  between  stability  and  ORE  is  illustrated  by 
comparison  of  the  ORE  for  NAVSTAR  1  at  0640  on  9  Hay  (fig. 
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T.t*)  «ith  the  stability  plot  on  figure  7.1  for  the  same  time 
period.  On  figure  7.2  the  stability  plot  variation  at  07U5 
on  10  May  remains  unexplained  but  may  have  been  due  to  a 
receiver  malfunction. 

Results  for  tests  on  11  May  (fig.  7.6^  7.7,  and  7.8), 
exhibited  less  stability  than  expected.  Also,  system  2  was 
less  stable  than  system  1.  The  largest  difference  in 
reading  between  the  two  receivers  occurred  at  0750  (fig. 
7.7).  This  difference  may  have  been  caused  by  the  location 
of  the  antenna  of  system  2  under  the  antenna  of  system  1. 
At  approximately  0800,  the  satellite  constellation  elevation 
angles  were  the  largest,  between  56  to  80°  above  the 
horizon  (fig-  7.9) ,  and  antenna  interference  may  have  led  to 
the  large  differences.  Neither  the  satellite  constellation 
azimuth  angles  (fig.  7.  10)  nor  the  user  range  errors  (fig. 
7.11)  appeared  to  have  any  effect  on  the  stability. 

Comparison  of  the  stationary  GPS  receiver  positions  with 
the  position  obtained  with  the  geoceiver  is  depicted  in 
figure  7. 12  for  the  stability  and  static  tests.  The 
comparison  indicates  that  there  is  a  nearly  constant  offset 
of  the  two  positions  with  the  GPS  position  almost  due  south 
of  the  geoceiver  position.    Further  comparisons  between  the 
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stationary  GPS  receiver  position  ani  the  geoceiver  position 
reinforce  this  conclusion.  In  figures  7.13,  7.1U,  and  7.15, 
the  mean  distances  between  the  two  positions  are  127,  173 
and  1U3  meters,  respectively.  The  higher  means  for  the  data 
shown  in  figures  7.1U  and  7.15  are  due  to  the  ORE  (fig.  7.16 
and  7,17)  for  those  dates  and  the  NA7STAR  upload  frequency. 
If  the  data  from  13  May  is  deleted  because  of  the  high  ORE 
for  that  date  (fig.  7.16),  and  the  median  distance  used  for 
15  Hay  is  taken  to  be  a  representative  value  (130  meters  on 
fig.  7. 15) ,  The  difference  is  approximately  129  meters 
Between  the  two  positions.  The  azimuth  was  also  nearly 
constant  with  the  observed  median  values  of  181°  ,  207°  and 
178  "  from  south  as  shown  in  figures  7.18,  7.19,  and  7.20 
respectively.  Neglecting  the  data  on  13  May  because  of  QRE 
values  (fig.    7.16),   the  remaining   two  values  are  nearly 

o 

identical  and  have  a  mean  of  179.5  from  south.  Thus  the 
data  indicates  a  nearly  constant  offset  between  the  position 
obtained  with  the  geoceiver  and  the  position  obtained  with 
GPS.  This  difference  is  about  130  meters  and  has  an  almost 
north-south  orientation. 
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Figure  7.13, 
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Figure  7.14. 
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Figure  7.15, 
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Figure     7.19, 
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Figure  7.20. 
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B.   E7AL0ATION  OF  N0NDI7FERENT lAL  GPS 

j      Evaluation  of  the   non-differential   GPS  positions   was 

I   done  on  a  set  by  set  basis.  Figures  7.21  through  7.31  depict 

1  the  distance  in  meters   between  the  theodolite-computed  ship 

position  and  the  shipboard  GPS  receiver  position.   Sets  1,  2 

and  3   have  the  smallest   distance  difference,   as   shown  in 

figures  7.21,  7,22^   and  7,23,   with  mean  values  between  118 

and  133  meters  and  standard  deviations  ranging  from  7.U5  to 

16  meters.   These  smaller  mean  distances  are  a  result  of  the 

low  user  range  errors  for  12   flay  as  shown  on  figure  7.32. 

The  highest  mean  distances  occurred  during  sets  U,   9,   10, 

12,   and   13  when  the  ORE  values  were  high  and  the  NA7STAR 

uploads  occurred  (fig.    7.16   and  7.17).    The   remaining 

distances  vary  between  1U4  and  278  meters. 

The  results  during  the  best  satellite  observation 
periods  are  depicted  in  figures  7.33,  7.34,  and  7.35,  where 
the  non-differential  GPS  ship  position  was  compared  to  the 
theodolite-computed  ship  position,  the  distance  between  the 
two  being  plotted  as  a  histogram.  Mean  distances  on  figures 
7.33  and  7.35  are  representative  of  the  five  satellites 
observed  and  the  improved  accuracy  of  the  system,  with  mean 
distances  of  134  and  149  meters  and   standard  deviations  of 
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Figure  7.21 . 
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Figure     7^22. 
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Figure  7.24. 
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Figure  7.25. 
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Figure     7.26 


♦*♦♦♦♦♦ 
1^  ill  in  11^  III  111  mf 

•   •••••• 


—Z— 

UJ<UJ 

o— o 
zoz 

XXX 


Z     aiiuuJUJiw 


3      ^^-Z^^ 

(XZOOOOOX 

'''ZO^/^OlnOX 

ox  •  •  •  •  ttS 


►-000000 

Z4.*    I     I    ♦♦ 

'  WW  f"  'in"  Miiin 
X  f^  in;  C^  ^  -O -^ 


vu 

O          i/K/1 

i/)»i 

CC               LUI/I 

uj       za— <M 

X        It— «■« 

O         ujac^-»- 

— <^-^«3uJUJ 

xxx</»5£<a(a 

2 

O 

M« 

►- 

<NJ.'<Mi.«i«^ 

OOOCDOO 

>/iZ 

**!♦♦* 

oo 

H  Mil  lint  ill  JUi 

a.>-i 

,r<of~<7'-oo 

>— 

^r- j-casrj 

a  — 

r-o-J'asiNOO 

— v> 

fMU^om'0*n 

xo 

fMO-*<^a<» 

>/lQ. 

00~»-UM»>«D 

IMOL 
t-X. 


UJt/1 
UiO. 

LUQ 


aloe 

ujO 

XO 

z 

Z3 

OJ 


zo 
-tz 
t-< 


a 

OXuj     ui 

azui>o    cc 

UJ—     u.ZOut 
of  OCOOI'IZQ 

v/>>i/loXflCX 


pgrMrsirMrsirsinj 
OOOOOOO 
-%  ♦  4-  ♦  •*  ♦  ♦ 

jUlULUIUUItUliJ 

»-/ r*  o  o  v/«  f^  Hi  »n 
Z '*^  ^ -f  t^tM  **^'^ 
jj   •••»••• 


3  >-  ZZ 

^  ZZOujUi 

-«      Z<<ZXX 

O  *£       'tujiij-4 

O-  i-z  — r'TaCTT 

UJ  .iJUJUJOC"**"Ui*t 

</<  OZXH-ZXOX 


164 


Figure  7.27, 
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Figure  7.28, 
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Figure  7.29, 
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Figure  7.30. 
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Figure  7.31  . 
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Figure  7.33. 
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Figure  7.34. 
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Fiqure   7.35. 
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10  and  52  meters.  The  large  standard  deviation  for  15  May 
is  apparent  in  the  bimodal  empirical  frequency  curve  (fig. 
7.35)  and  is  caused  by  the  NAVSTAR  upload  (fig.  7.17).  in 
this  case,  the  GPS  receiver  did  not  obtain  the  new  ephemeris 
until  halfway  through  the  time  period  of  the  observations. 
The  data  for  13  May  (fig.  7.3U)  reflects  frequent  NAVSTAR  1 
uploads  and  the  resulting  high  URE  values  (fig.  7.16)  for 
that  date. 

C.   DIFFERENTIAL  GPS  EVALaATION  WITH  STATIC  CORRECTION 

The  results  of  the  evaluation  of  the  data  of  12  May 
(sets  1,  2,  and  3),  computed  using  the  differential  GPS 
positions  and  the  static  corrections,  are  shown  in  figures 
7.36,  7.37,  and  7.38.  Mean  difference  between  positions  for 
the  theodolite  ship  position  and  the  GPS  position  are  29.1, 
27,4,  and  64.8  meters,  with  standard  deviations  of  7.7,  7.5, 
and  22.7  meters.  These  distances  are  larger  than  those 
determined  for  the  theodolite  ship  position  and  the 
differential  GPS  position  (fig.  7.39,  7.40  and  7.41).  This 
higher  mean  distance  indicates  that  the  static  correction, 
as  applied,  results  in  a  larger  error  in  the  differential 
GPS  position.  Therefore,  further  evaluation  of  differential 
GPS  with  static  corrections  were  not  conducted. 
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Figure  7.36. 
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Figure  7.40 
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Figure   / .41 . 
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D.      EVALUATION   OF    DIFFERENTIAL    GPS   WITHOUT    STATIC 
COHHECTIONS 

Evaluation  of  differential  GPS  without  the  static 
corrections  were  computed  by  sets  for  the  five  satellite 
constellations  for  each  day  of  testing.  Set  1  (fig.  7.39) 
shows  a  mean  distance  of  10.1  meters  between  the  two 
positions  and  a  standard  deviation  of  5.0  meters.  The 
median  azimuth  between  the  two  was  257°  from  south  (fig. 
7.42),  which  also  was  the  approximate  course  of  the  vessel 
during  the  test.  These  values  indicate  that  the  error  is 
caused  by  the  bias,  or  alias,  in  GPS  position  due  to  the 
Manpack  2U-second  fixing  interval  and  positioning  algorithm. 
For  this  set,  the  GPS  position  was  obtained  after  the 
theodolite  position. 

During  set  2  the  mean  distance  difference  between  the 
two  positions  was  11.7  meters,  and  the  standard  deviation 
was  6.8  meters  (fig.  7.40);  the  median  azimuth  was  71 
(fig.  7.43)  .  These  values  were  slightly  greater  than  those 
of  set  1  due  to  the  increased  URE  for  NAVSTAR  1  during  the 
set  (see  0700  on  fig.  7.32)  Additionally,  the  azimuth  was 
nearly  180°  from  that  computed  in  set  1,  with  the  vessel  on 
the  same  track  as  set  1.   Unlike  the  error  obtained  from  the 
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Figure  7.42. 
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Figure  7.43. 
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Manpack  bias  on  set  1,  the  error  in  set  2  was  due  to  a 
similar  time  delay  bias  because  the  theodolite  position  was 
obtained  after  the  GPS  position. 

Set  3  values  on  figures  7.41  and  7.44  are  similar  to 
those  of  set  2,  except  the  mean  distance  values  were 
greater — 44,7  meters  for  set  3  versus  11.7  meters  for  set  2. 
Mean  distance  increases  for  set  3  were  due  to  the  upload  of 
NA7STAR  1  (fig.  7.32)  and  the  availability  of  only  a  four 
satellite  constellation.  Set  2  had  a  five  satellite 
constellation  available. 

For  set  4,  figures  7.45  and  7.46  show  the  distances  and 
azimuths  for  the  set.  These  histograms  fail  to  illustrate 
the  constantly  increasing  distances  between  the  two 
positions.  These  increasing  distances  result  from  the 
uploads  of  NAVSTAR's  5  and  6  and  the  large  ORE  for  MAVSTAR  1 
(fig.  7. 16) .  Onboard  the  vessel  the  new  ephemaris  data 
uploaded  for  NAVSTAR 's  5  and  6  was  not  obtained  until  late 
in  the  set.  Also  near  the  end  of  the  set,  NAVSTAR  1  was 
uploaded  with  the  result  that  the  last  several  distances 
were  very  large,  as  shown  in  figure  7,16. 

In  set  5,  a  very  high  QRE  value  for  NAVSTAR  1 
contributed   to  the   extremely   large   position   difference 
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Figure  7.44 
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Fig-ire   7.4F 
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Figure  7.46. 
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between  that  determined  by  the  theodolites  and  that  for  the 
differential  GPS  position  (fig.  7.U7).  The  azimuth  remained 
nearly   constant  between   the  two   positions   with  a   median 

o 

azimuth  of  43  (fig*  7.U8).  These  values  for  set  5  show 
that  the  large  ORE  for  NAVSTAR  1  was  the  most  important 
contributing  error  source  (fig.  7,16).  The  large  ORE  may 
have  been  caused  by  activation  of  the  hydrazine 
stabilization  jets  to  dampen  momentum  on  NA7STAR  1. 

For  set  6,  the  mean  distance  difference  between  the  two 
positions  was  60.0  meters  and  a  median  azimuth  of  93  (fig. 
7.49  and  7.50).  These  values  indicate  that  the  upload  of 
NAVSTAR  1  near  the  beginning  of  the  set  (0715  on  fig.  7.16) 
was  received  and  utilized  in  the  Hanpack  position  algorithm, 
but  the  higher  distance  between  the  positions  also  indicates 
a  large  ORE  (also  shown  on  fig.  7.16). 

Data   from   set   7   were    not  evaluated  because   the 

stationary  ManpacJc  receiver  did  not  record  the  GPS  position 

for  each  fix   due  to  radio   failure.     The  mean   distance 

,  between  the  two  positions  for  set   8  was  121  meters   with  a 

I  mean  azimuth   of  63  °     (fig.    7.51   and  7.52).    Both   the 

„  distances  and   azimuths  show  that  the   URE  for  NA7STAR   1  at 

I 

0530   (fig.   7.17)    was  increasing   in   magnitude  and   thus 
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Figure  7.47. 
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Figure  7.48. 
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Figure  7.49. 


O  I 
Ol 


ri  III  III  m  111  111  iiji 

C(*^-"U^ol'^0 


LU4UJ 

zoz 

IXX 


2>      ZZZ2Z 

xa:<<<'<«ix 

—3333333 
aCXOOOOOX 

or X 


►-000000 
Z  ♦  ♦  ♦  ♦  ♦  ♦ 

UjUJLUUjUJUJUi 

C  — — O^»^OP• 
X<3Dtf^^rr*^0^U^ 

o    z 

>-       UJ 

ULI  VO»M 

IX    3e       UJU1 

«      UJ  ZO-<*>« 

U       X  !•-<< 

— <«^  >»■■)£  3  luui 

T*-    X7X</*>£aica 

03 


ax 


ir:?:  u-irrr  3:1:^3:3:1:  s:3:rxrjmr-z3:3:Trxi:T-mrmj:3:rrrmTi:3:x 


000000 
a.  ♦♦!♦♦♦ 

Z  UJtUUUUJUIUJ 

»/ifc«         p^f^-r^OW 

<mO         (r\0*^<©(*^'0 
I  Q.  ....«• 

1/1 

ZQ. 

luo  O 

UJ  UJ     CO     < 

oxuj     Ui 

aziij>o    cc 

S^.  UJ»«     u.ZO(/l 

Z  iXqlOuj<ZO 

uouj  a.»th-Ouj<*^ 

cz.(x  i/>  >i/ivj  zoe  X 


2_l 


i    "-^ 


'♦«*«♦  ^ 


3  ooooooo 

I  *JJ  *  *  *  ^  «  «  « 

X  UJIUUJUJUJUiUJ 

■yi  uj««**»«« 


^  zz 

^  ZZOUiUJ 


19! 


Figure  7.50. 
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Figure  7.51 . 
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Figure  7.52 
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resulting  in  a  greater  positional  difference.  Another 
factor  was  the  bias  introduced  by  the  Manpack  2U-second 
reception  period.  This  is  evident  in  the  azimuth  for  the 
theodolite  ship  position  obtained  from  the  GPS  shipboard 
receiver. 

The  data  for  sets  9  and  10  were  combined  for  evaluation 
and  the  results  were  a  mean  distance  between  positions  of 

o 

241  meters  and  an  azimuth  of  267  (fig.  7.53  and  7.5U). 
Values  of  that  magnitude  show  that  the  ORE  for  NAVSTAR  1  was 
the  error  source.  Data  were  sampled  from  0600  to  0640  (fig. 
7.17).  Thus  the  correlation  between  the  high  NA7STAR  1  ORE 
and  the  high  mean  distance  difference  between  positions. 
Again,  receiver  bias  also  introduced  an  error  into  the 
evaluation. 

Analysis  of  set  11  shows  a  bimodal  distribution  in  both 
the  two  position  differences  and  azimuths  (fig.  7.55  and 
7.56).  The  peaks  on  the  distance  histogram  are 
approximately  133  and  160  meters;  the  peaks  on  the  azimuth 
histogram,  approximately  220  °  and  236° .  Both  peaks  are 
caused  by  the  upload  of  NA7STAR  3  halfway  through  the  set 
and  the  increasing  magnitude  of  the  MAVSTAR  1  ORE  just  after 
the  upload  as  shown  by  figure  7.17.  The  total  positional 
difference  for  set  11  was  approximately  154  meters. 
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Figure  7.53. 
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Figure  7.54. 
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Figure  7.55. 
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Figure  7.56. 
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Evaluation  of  the  combined  sets  12  and  13  (fig.  7.57  and 
7.58)  result  in  a  mean  distance  between  the  two  positions  of 
172  meters  and  a  mean  azimuth  of  262° .  The  mean  distance 
was  large  because  of  the  large  magnitude  of  the  NA7STAR  1 
ORE  (fig.  7.17).  A  contributing  source  of  error  was  the 
bias  introduced  by  the  Manpack  receiver. 

Overall,  the  data  using  a  five-satellite  constellation 
were  as  good  or  better  than  the  data  using  only  four 
satellites.  The  results  for  12  May  using  five  satellites 
(fig.  7.59)  show  a  mean  distance  of  11.4  meters  and  a  median 
of  10.6  meters.  This  compares  favorably  with  the  results 
from  sets  1  and  2  (fig.  7.39  and  7.40),  which  had  mean 
distances  of  10.1  and  11.7  meters,  respectively.  The  five 
satellite  data  from  13  Say  showed  a  very  large  variance  due 
to  the  large  NA7STAR  1  ORE.  Therefore,  data  were  not 
compared  to  remaining  13  May  data.  Data  on  15  May  for  the 
five  satellite  constellation  had  a  bimodal  distribution 
because  NAVSTAR  1  was  uploaded.  Once  the  uploading 
occurred,  the  peak  distance  was  137  meters  (fig.  7.50)  which 
compares  favorably  with  121,  2U1,  154,  and  172  meters  for 
sets  a,  9  and  10,  11,  12  and  13.  Thus  with  a  five-satellite 
constellation,  positioning  results  were  as  good  or  better 
than  the  four-satellite  constellation. 
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Figure  7.57 
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Figure  7.58. 
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Figure  7.59. 
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Figure  7.60. 
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Figures  7.61  to  7.74  show  the  diffarence  by  fix  between 
the  latitude  and  longitude  plots  of  the  theodolite  ship 
position,  the  differential  GPS  position,  and  the  uncorrected 
GPS  position.  During  most  of  the  testing,  the  longitude  of 
both  methods  of  GPS  positioning  was  in  closer  agreement  with 
the  theodolite  ship  position  than  was  the  latitude.  Shown 
on  figure  7.9  and  7.10  are  the  satellite  elevation  and 
azimuth  data  for  the  testing  period.  This  data  indicates 
that  the  satellite  constellation  had  a  larger  effect  on 
latitude  than  longitude  because  the  predominant  ranging  was 
computed  from  one  quadrant. 

The  large  differences  between  the  GPS  position  relative 
to  the  theodolite  ship  position  were  easily  identified 
during  the  shipboard  testing  period.  Latitude  and  longitude 
plots  (fig.  7.65,  7.66,  7.69  through  7.74)  for  sets  4,  5,  9, 
10,  and  12  show  the  deviations  of  the  GPS  receiver  positions 
from  the  predicted  positions  that  were  obvious  to  the  GPS 
operators  during  the  shipboard  testing.  For  sets  4  and  5, 
the  stationary  GPS  receiver  positions  began  to  deviate 
greatly  from  the  geoceiver  position,  resulting  in  a  large 
differential  correction  to  the  GPS  ship  position.  These 
large  differential   corrections  caused  the  differential  GPS 
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position  shown  on  figures  7.6  5  and  7,66  to  deviate  a  large 
amount  from  the  theodolite  position  and  the  uncorrected  GPS 
position. 
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7III.   CONCLDSIONS 

The  objective  of  the  study  was  to  determine  whether  the 
application  of  GPS  in  a  translocated  or  differential  mode 
provides  sufficient  position  accuracy  for  near-shore 
hydrographic  surveys.  Several  assumptions  were  made  prior 
to  the  test:  (1)  two  GPS  receivers  required  for  the 
differential  mode  performed  the  computations  similarly,  so 
that  in  the  worst  case  a  constant  bias  would  be  displayed  if 
both  receivers  were  colocated;  (2)  the  process  of  updating 
the  ephemeris  information  would  not  contribute  a  significant 
positional  error  to  GPS  in  the  differential  mode;  (3) 
reasonable  accuracy  could  be  obtained  by  using  the  output 
from  the  receiver  control  display  unit,  that  is, 
sophisticated  processing  methods  were  not  required;  and  (U) 
the  Manpack  receiver  results  would  not  differ  when  used  in 
stationary  or  low-dynamic  environments  (shipboard  at  U 
knots)  . 

None  of  these  assumptions  proved  to  be  valid.  It  was 
originally  expected  that  the  two  receivers  would  produce 
similar  positions  if  both  sets  were  operating  at  the  same 
time  over  the  same  location.  One  limitation  of  these 
receivers,  however,  is  that  current  ephemeris  information  is 
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not  received  unless  the  receiver  has  been  activated  prior  to 
the  update.  [Re£.  66]  For  a  period  of  time  after  the 
receiver  is  activated,  the  portion  of  the  satellite  message 
that  contains  the  updated  ephemeris  is  not  acknowledged. 
Therefore,  position  degrades  rapidly  when  a  NA7STAR  is 
updated  and  the  current  ephemeris  information  is  not 
utilized  by  the  receiver.  For  this  study,  unfortunately, 
the  initial  times  of  activation  of  the  two  receivers  were 
not  recorded.  Hence,  the  effect  on  position  caused  by  this 
receiver  parameter  was  not  defined. 

The  ORE  plots  illustrate  that  the  ephemeris  updates  did 
contribute  significantly  to  the  range  error  and, 
consequently,  the  resulting  position  error.  The  position 
determined  by  the  Manpaclc  receivers  was  not  affected  as 
seriously  by  the  clock  error  in  NAVSTAR  1  when  a 
five-satellite  constellation  was  used.  The  inaccuracy  of 
pseudo- ranges  received  from  this  satellite  required  frequent 
updates  from  Vandenburg,  especially  on  13  May.  A  receiver 
that  had  the  option  of  totally  eliminating  the  erroneous 
pseudo-range  would  have  given  the  best  results  for  the 
study.  Since  the  receivers  operated  independently,  and  the 
ephemeris  information  was  not  concurrent,  the  differential 
mode  was  effective  for  only  a  few  cases. 
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A  simple  correction  was  applied  to  the  values  output  on 
the  control  display  unit  of  the  Manpack  receiver.  This  did 
not  permit  application  of  tropospheric  or  refraction  error 
corrections,  other  than  that  applied  by  the  receiver  during 
computation  of  the  position,  &  more  sophisticated 
processing  method,  resulting  in  greater  accuracy,  would 
require  acquisition  of  the  pseudo- ranges  and  subsequent 
postprocessing  to  determine  the  receiver  positions. 

The  best  results  during  the  onboard  portion  of  the  study 
were  obtained  when  the  ship  was  allowed  to  drift  at 
approximately  1  knot  on  12  May. 

In  conclusion,  use  of  two  Phase  I  Manpack  receivers  in 
the  differential  mode  meets  the  required  accuracy  of  10 
meters  only  if  the  ship  speed  is  less  than  1  knot  and  other 
parameters  such  as  ephemeris  information  and  the  number  of 
satellites  used  for  the  position  computation  introduce  a 
minimum  amount  of  error.  Further  testing  using  differential 
GPS  should  incorporate  more  sophisticated  processing  methods 
and  more  advanced  receivers  to  achieve  the  survey  accuracy 
required  for  nearshore  hydrographic  surveys. 
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APPENDIX  A^  DESCRIPTION  OF  ERROR  THEORY 
From  applicatioa  of  error  theory  in  the  evaluation  of 
positional  information,  it  is  possible  to  establish  a 
meaningful  accuracy  statement  subject  to  uniform 
interpretation.  To  provide  a  logical  and  acceptable  basis 
FOR  COMPOTATION  AND  COMPARISON,  POSITIONAL  ERRORS  ARE 
assumed  to  follow  a  normal  distribution.  The  assumption  is 
valid  because  positional  error  components  generally  follow  a 
normal  distribution  pattern  when  sufficient  data  is 
available. 

The  statistical  treatment  of  errors  utilizes  measurable 
errors  obtained  in  the  sources  of  positioning  information. 
Analysis  of  the  linear  components  provides  a  two-dimensional 
expression  of  the  accuracy  of  the  positioning  system  in  the 
form  of  an  error  ellipse.  The  use  of  an  error  ellipse  is 
complicated  by  the  problem  of  axes  orientation  and 
propagation  of  elliptical  errors.  Therefore,  the  ellipse  is 
commonly  replaced  by  circular  form  which  is  easier  to  use 
and  understand.  The  linear  standard  errors  are  combined  and 
converted  to  a  circular  distribution  and  the  final  accuracy 
statement  is  expressed  in  terms  of  circular  errors. 
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The  circular  form  assumes  the  circular  probability 
distribution  function.  This  function  expresses  the 
probability  that  the  radial  error  will  be  equal  to  or  less 
than  the  radius  of  the  probability  circle.  The  dispersion 
of  errors  within  the  distribution  is  measured  in  terms  of 
precision  indexes,  each  of  which  represents  the  error  which 
is  unliJcely  to  be  exceeded  for  a  given  probability.  The 
preferred  circular  precision  indexes  are  the  circular 
standard  error  {o  )  ,  the  circular  probable  error  (CPE  or 
CEP),  the  circular  map  accuracy  standard  (CMAS) ,  and  the 
circular  near-certainty  error,  or  3.5  sigma  (3.5  a  ).  For 
each  precision  index,  the  given  probability  is  39.35%,  50%, 
90%,  and  99.78%,  respectively. 

Each  of  the  precision  indexes  can  be  converted  to  one  of 
the  others  by  using  the  following  table: 

Circular  Error  Conversion  Factors 


39.35% 

50% 

90% 

99.78% 

39.35% 

1.0000 

1.1774 

2.1460 

3.5000 

50% 

0.8493 

1.0000 

1.32  27 

2.9726 

90% 

0.466  0 

0.5486 

1,00  00 

1.6309 

99.78% 

0.2857 

0.3364 

0.6131 

1  .0000 
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The  rapid  approximation  for  deteraining  the  circular 
standard  error  is 

CSE  =  0.5000  (ax  ♦  ay) 
where  f^  x  and  ^y  are  linear  standard  errors  obtained  from  the 
sources  of  positioning  information,  and  therefore  are  the 
axes  of  the  error  ellipse.  For  a  truly  circular 
distribution,  where '^  c  is  identical  to  the  circular  standard 
error,  a  x  and  a  y  are  equal  and  the  angle  of  intersection 
between  the  two  is  90°.  In  all  other  cases,  a  normal 
circular  error  distribution  is  substituted  for  the 
elliptical  distribution.  The  substitution  is  satisfactory 
for  an  analysis  within  specifiedcmin/omax  ratios,  where  a  min 
is  the  minor  or  smaller  linear  standard  error  of  the  two. 
Distortion  in  the  error  distribution  berween  elliptical  and 
circular  forms  at  lower  a min/ a  max  ratios  limits  the  use  of 
the  circular  concept  to  ratios  greater  than  0.2  . 

The  CMAS  is  the  precision  index  used  as  the  O.S. 
National  Map  Accuracy  Standard.  It  is  interpreted  as 
limiting  the  size  of  error  which  90%  of  the  positions  will 
not  exceed.  CMAS  is  computed  from  circular  standard  error 
as  follows: 

CMAS  =  2.  1U60  a  c 
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The  two-dimensional  error  distributions  for  positions 
determined  by  the  ARTEMIS  system  and  the  theodolite  network 
are  described  by  CMAS.  The  accuracy  of  the  specific 
positions  is  expressed  by  a  statement  of  probability  and 
error  magnitude.  The  accuracy  statement  does  not  mean  that 
the  error  in  position  is  exactly  the  value  shown  in  figure 
5.1,  rather  it  expresses  the  probability  that  the  error  in 
position  will  not  be  larger  than  the  error  given. 
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71.  Commanding  Officer 
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Naval  Ocean  Research  and  Development  Activity 
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National  Ocean   Survey    (OA/C) 

National  Oceanic   and   Atmospheric   Administration 
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National  Oceanic   and    Atmospheric   Administration 
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National  Oceanic  and  Atmospheric  Administration 
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88.  Director  1 
Instituto    Hidrograf ico 

Rua   Das  Trinas,    49 
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